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Abstract 

This paper describes the cutting sequences of geodesic flow on the modular surface ^/PS'L(2, Z) 
with respect to the standard fundamental domain !F = {z = x + iy : —^<x<^ and \z\ > 1} 
of P5L(2,Z). The cutting sequence for a vertical geodesic {9 + it : t > 0} is related to a one- 
dimensional continued fraction expansion for 6, called the one-dimensional Minkowski geodesic 
continued fraction (MGCF) expansion, which is associated to a parametrized family of reduced 
bases of a family of 2-dimensional lattices. The set of cutting sequences for all geodesies forms 
a two-sided shift in a symbol space {L,R,J} which has the same set of forbidden blocks as 
for vertical geodesies. We show that this shift is not a sofic shift, and that it characterizes 
the fundamental domain up to an isometry of the hyperbolic plane f). We give conversion 
methods between the cutting sequence for the vertical geodesic {9 + it : t > 0}, the MGCF 
expansion of 9 and the additive ordinary continued fraction (ACF) expansion of 9. We show 
that the cutting sequence and MGCF expansions can each be computed from the other by a 
finite automaton, and the ACF expansion of 9 can be computed from the cutting sequence for 
the vertical geodesic 9 + it hy a. finite automaton. However, the cutting sequence for a vertical 
geodesic cannot be computed from the ACF expansion by any finite automaton, but there is 
an algorithm to compute its first i symbols when given as input the first 0{i) symbols of the 
ACF expansion, which takes time 0(^^) and space 0{i). 
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1. Introduction 

This paper describes the symbohc dynamics of cutting sequences of geodesies for the stan- 
dard fundamental domain = {z = x + iy : — ^ < x < ^, \z\ > l}of the modular group 
PSL{2, Z) acting on the upper half plane = {z G C : Im{z) > 0}. 

We study in particular the cutting sequences of vertical geodesies {6 + it : t > 0} for G M. 
These particular geodesies are related to a continued fraction expansion introduced in 1850 by 



case the family associated to 9 is 



Hermite [^0[ in terms of quadratic forms, and studied by Humbert [22, 23]. Our motivation 
for studying them was that they appear in the one-dimensional case of a multidimensional 
continued fraction introduced in Lagarias (3^. This expansion, called the Minkowski geodesic 
continued fraction expansion (MGCF expansion), is based on following a parametrized family 
of lattice bases in GL{d + 1, 'L)\GL[d + 1, M) as the parameter t varies; in the one-dimensional 

^ . The MGCF expansion has the merit that it finds 

— a t 

good Diophantine approximations in all dimensions. However, the one-dimensional case of this 
expansion does not coincide with the ordinary continued fraction expansion. We show in §3.4 
that the one-dimensional MGCF expansion of is essentially equivalent to the cutting sequence 
expansion of the vertical geodesic {d + it : Q <t < oo}. We use this connection to determine the 
precise relation of the MGCF expansion of 9 to the additive continued fraction expansion of 9; 
the additive continued fraction is the variant of the ordinary continued fraction which includes 
all intermediate conver gents. 

A second motivation for studying cutting sequences on ^/PS'L(2,Z) arises from their use 
as symbolic codings of geodesies in the study of geodesic flow on constant negative curvature 
surfaces of finite volume; we give more background on this at the end of the introduction. Cut- 
ting sequence encodings are of special interest because they apparently encode more information 
about the geodesic fiow on a Riemann surface than other symbolic encodings of geodesies, which 
have a simple description (shift of finite type) but only retain topological and not conformal 
information about the Riemann surface. Adler and Flatto j^, §10] raised the question whether 
cutting sequence encodings preserve conformal information about the Riemann surface and also 
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retain information about the lines on the Riemann surface used to define the cuts, i.e. the shape 
of the fundamental domain in the universal cover. Consider any hyperbolic polygon V in Sj 
which is a fundamental domain for a properly discontinuous group T acting on Sj, such that Sj/F 
has finite volume, and suppose that V has {Vl sides. Let S-p denote the closed subshift of the 
shift on \V\ letters generated by the cutting sequences for V for geodesic flow on general Sj/T, 
as defined in §2.2. Adler and Flatto |^, p. 3OO]0 ask: does S-p determine V up to conformal 
isometry and T up to hyperbolic conjugacy? We show here that this is the case for the standard 
fundamental domain of PS'L(2,Z). 

Our main results on cutting sequence encodings are as follows. 

(1) The set of cutting sequence expansions for (irrational) vertical geodesies 11^ is characterized 

in terms of forbidden blocks (Theorem 6.1). A generating set of forbidden blocks is 
enumerated. (Theorems and ^^■)- The number of minimal forbidden blocks of length 
at most k grows exponentially in k. (Theorem 6.3). 

(2) The cutting sequence expansions for all geodesies comprise (essentially) a closed two-sided 

subshift SjF of the full shift on three symbols {L, R, J}. The set of forbidden blocks of Sjf 
is the same as that for 11^. (Theorem 7T). Each symbol sequence in corresponds to 
a unique oriented geodesic on Sj/PSL{2,Z,), with the exception of the two sequences L°° 
and R°°. In the converse direction, each oriented geodesic gives rise to a finite number 
of shift-equivalence classes of symbol sequences in Sjp. This number is at most eight if it 
is not a periodic geodesic. (Theorem 7.2). The shift T,jr is not a sofic system. (Theorem 
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(3) The shift Ejf characterizes J- up to an isometry of i^. If "P is a hyperbolic polygon of 
finite area (possibly with some ideal vertices) which is a fundamental domain of a discrete 
subgroup r of PSL{2,M.) and if the subshift S-p is isomorphic to Sjf by a permutation 
of symbols, then there is a hyperbolic isometry g G PSL(2,M.) such that = gV and 
PSL{2,Z) = gVg-^. (Theorem O) . 



fractional transformation z — > Nz with N 



In obtaining result (1) we show in §5 that vertical geodesies have special features (not shared 
by general geodesies) that facilitate characterizing forbidden blocks: any vertical geodesic that 
hits a PSL{2, Z)-translate of a corner of the fundamental domain must have rational 9, and 
each such geodesic hits at most one such corner, with the exception of those 9 = ^(mod 1), 
which hits two corners. Furthermore, the continued fraction expansion of any rational 9 such 
that {9 + it : t > 0} hits a corner has a symmetry which can be described in terms of the linear 

1 2 

2 1 

(3) we define cutting sequence expansions for geodesies that hit a corner of the fundamental 
domain of S)/PSL{2,'Z) to be limits of general-position geodesies that hit no corner. This 
procedure assigns infinitely many different cutting sequences to certain periodic geodesies. We 
interpret result (3) as showing that the shift Sjf encodes conformal information about J^. 

As a preliminary to obtaining the above results, in §2-§5 we determine precise relations of the 
cutting sequence of the geodesic 9 + it to the Minkowski geodesic continued fraction expansion 
of 9 and to the additive version of the ordinary continued fraction expansion of 9, which we call 



see Theorem 5.1. In formulating result 



^They raise the question for the {8g — 4)-sided fundamental polygons discussed in j^, for genus g > 2. 
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the ACF expansion. (The ACF expansion, which is described in §3.1, includes all intermediate 
convergents.) These give algorithms for converting between the symbolic expansions, and obtain 
results on how hard it is computationally to convert between these different symbolic expansions, 
as follows. 

(1) The cutting sequence expansion of {9 + it : t > 0} and the Minkowski geometric continued 
fraction expansion of 9 can each be computed from the other by a finite automaton. 



(Theorem 3^). The additive ordinary continued fraction expansion of 9 can be computed 



from either of these expansions by a finite automaton. (Theorems 4.1 and 4.2) 



(2) The Minkowski geodesic continued fraction expansion of 9 cannot be computed from the 

additive ordinary continued fraction expansion of by a finite automaton. (Theorem |6.4| ). 

(3) The Minkowski geodesic continued fraction expansion of 9 can be computed from the 
additive ordinary continued fraction expansion of 9 in quadratic time and linear space on 
a random access machine, i.e. the first I MGCF digits can be computed in time 0{i'^) 
and space 0{i), from 0{£) symbols of the ACF expansion. (Theorem |5T 



The notion of finite automaton (finite state machine) is described in §3.5. The key feature 
of a finite automaton is that it has a fixed finite amount of memory. Result (1) requires that 
given the successive symbols in the input expansion there is a an absolute bound B so that 
after scanning B consecutive input symbols at least one new output symbol is determined. This 
property may be called bounded look-ahead. Result (2) comes from the fact that the amount of 
look-ahead needed to compute one MGCF symbol from the ordinary continued fraction expan- 



sion can be unbounded. Result (3) is based on a conversion algorithm given in Theorem 5.1 
To compute a given output symbol it potentially requires remembering the entire input string 
of symbols to that point, which is the source of the linear space requirement. This computation 
can be carried out in the given time and space in the random access machine(RAM) model for 
computation, which allows unbounded storage. It can also be carried out on a one-tape Turing 
machine using 0{i) space with a time bound polynomial in £ on a one-tape Turing machine. 
For the random access machine (RAM) and Turing machine models of computation, see Aho, 
Hopcroft and Ullman Q. Results (l)-(3) delineate the computational relations between the or- 
dinary continued fraction expansion, the MGCF expansion and the cutting sequence expansion. 
Result (2), which is the significant result here, shows that cutting sequences encode information 
more concisely than the additive ordinary continued fraction expansion. 

To conclude this introduction, and to put the results above in a more general context, we 
recall background on symbolic codings of geodesic fiow on a compact Riemann surface with a 
finite number of punctures. If F is a finitely generated discrete group of conformal isometries 
of the hyperbolic plane Sj, such that Sj/T has finite volume, then Sj/T is a compact Riemann 
surface^ minus a finite number of punctures, and geodesic fiow on Sj/T contains both topological 
and conformal information about this Riemann surface. Symbolic codings of geodesies were 
introduced by Hadamard |l^ in 1898 as a way of understanding the complicated motions 
on geodesies on such surfaces. E. Artin [0] showed in 1924 that there exist dense geodesies 



^There are some mild extra conditions needed on F for compactness, see Lehner [ p3| , pp. 203-205] for sufficient 
conditions. It suffices for T to have a fundamental region T that is a hyperbolic polygon with a finite number of 
sides, with only parabolic vertices at infinity. 
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on the modular surface S)/PSL(2,'Z), using symbolic encodings with continued fractions, and 
in 1935 G. Hedlund |18| used Artin's coding to show that geodesic flow on this surface was 
ergodic. Some other symbolic encodings were "boundary expansions" by Nielsen [41| and 
cutting sequences with respect to a fundamental domain of T by Koebe |2?] and Morse p9[ . 
Cutting sequence expansions can be viewed as a generalization of the reduction theory of 



indefinite binary quadratic forms to arbitrary finitely-generated Fuchsian groups, see Katok [24|. 
More recently Bowen and Series @ and Series ||46| , |47|| used "modified boundary expansions ' 
to obtain a particularly simple symbolic expansion from which strong forms of ergodicity could 
be deduced; their codings are sofic systems. Adler and Flatto |2|, |^ used "rectilinear map" 
codings of geodesic flow on Sj/T to explain why certain specific maps of the interval, such as 
the continued fraction map and backwards continued fraction map, have invariant measures 
of a simple form. They showed that these invariant measures are inherited from the invariant 
measure on geodesic flow (Liouville measure) using a cross-section map followed by a factor 
map. The codings of Adler and Flatto Q have a particularly simple structure - they are shifts 
of finite type - and they preserve topological information about the Riemann surface, but they 
lose conformal information, for they are identical for all Riemann surfaces of the same genus 
g > 2. (See |^, Theorem 8.4 and Section 10].) In this respect it is of interest whether cutting 
sequence encodings preserve conformal information. Besides all these encodings, two other 
encodings of geodesic flow on the modular surface S^/PSL{2,Z) with interesting Diophantine 
approximation properties appear in Arnoux H] and Lagarias and Pollington |31]. 

With regard to this general framework, we note that the proofs in this paper heavily rely 
on specific facts about PSL{2,'Z). The methods used may extend in some form to congruence 
subgroups of PSL{2,Z), i.e. arithmetic Fuchsian groups, but do not appear likely to apply to 
general finitely generated Fuchsian groups F. 

For general terminology in symbolic dynamics we follow Adler and Flatto |p and Lind and 
Marcus |^^. From the viewpoint of symbolic dynamics, the associated shift Sjf studied in this 
paper is an example of a naturally occurring shift more complicated than the ones currently 
having a well-developed theory, cf . l33| . 



2. Cutting Sequences 

We describe cutting sequences for finite-sided hyperbolic polygons that are the fundamental 
domain of a given discrete subgroup F of PSL{2,M) that acts discontinuously on Sj, and then 
specialize to the standard fundamental domain of the modular group PSL{2,7j). General 
references for geodesic fiow and for cutting sequences are Adler and Flatto 0] and Series |49]. 



2.1. Cutting Sequence Shifts 

Let Sj denote the hyperbolic plane, represented as the upper half-plane ?) = {z = x + iy : 
y > 0} with hyperbolic line element ds'^ = ^{dx^ + dy'^) and volume An oriented geodesic 
7 = (01, 62) in ^ is uniquely specified by its two ideal endpoints 61,62 G MU {00}, with 61 / 62- 
It is oriented from its head 61 to its foot 62- A vertical geodesic (00, 9) is the vertical line 6 + it, 
with t decreasing from cxd to 0. Other geodesies are semicircles perpendicular to the real axis 
with endpoints 61 and 62, centered at and oriented to start at 61 and travel to 62. 
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Let r be a finitely generated discrete subgroup of PSL(2, M) which acts properly discontinu- 
ously on the upper half-plane Sj = {x+iy : y > 0} as a group of linear fractional transformations. 
We suppose that Sj/T has finite volume with respect to the hyperbolic volume on S^. Let V be 
any finite-sided convex hyperbolic polygon which is a fundamental domain for T, i.e. V is such 
that the set of F-translates {gV : g £ T} tessellates S), and no element g G T sends V to itself 
except the identity. We allow polygons V to have one or more ideal vertices ( "cusps" ) , and they 
may also have extra vertices, called elliptic vertices, which have an angle of vr and which fall 
in the middle of a side of V regarded as a geometric object, as discussed below. Fundamental 
domains exist for each such finitely generated F, cf. [4, Theorem 10.1.4]. Such polygons V 
necessarily have an even number of "sides," where by convention each elliptic vertex divides 
a side of the geometric object V into two "sides." If V is to be a fundamental domain for F, 
the "sides" of V must be identified under the action of F. Specifically, there is a pairing of the 
"sides" of V which assigns to each "side" s of V a "side" s' and an isometry g{s,s') G F such 
that g{s,s') maps s to s' . Furthermore (s')' = s and 

9{s',s) = {g{s,s'))~^, 

while if s = s' then g{s,s) is the identity on s. The pairing indicates how V tessellates Sj. 
For g = g{s,s'), the domain of V is adjacent to V in such a way that the "side" labeled s of 
gV coincides with the "side" labeled s' of V, while g{Int(V)) n Int{V) =0. If s = s' then 
necessarily g^ is the identity on S), and is a hyperbolic reflection. The neighborhood set of V 

Nr{V) ■={geT ■.g = g{s, s') for some side s of V} . (2.1) 

is a set of generators of the group F. It will comprise the set of symbols used in the cutting 
sequence encodings of geodesies described below. 

The number of elements of Ny{V) is generally equal to the number of "sides" of V, but can 
be less if there are two or more elements g{s,s') that are equal. For convex polygons V this 
situation can only arise when there are hyperbolic reflections g = g{s, s') with s ^ s' , in which 
case g{s, s') = g{s' , s). In such a case the "sides" s and s' must lie on a common geodesic of ^ 
and share a common endpoint,^ which is then a vertex of V having an angle of vr. We call any 
such vertex an elliptic vertex. These two "sides" s and s' are then assigned the same generator, 
so if we erase the elliptic vertex and glue them together, we get a single geometric side of 
the polygon V regarded as a geometric object. We therefore have: For a convex fundamental 
domain "P of a finitely generated group F, the number of elements of Nr{P) is equal to the 
number of geometric sides of V. In particular, the number of elements of Nr{P) may be odd. 

As an example, the standard fundamental domain of T = PS'L(2,Z), is geometrically a 
triangle with one ideal vertex, but as a fundamental domain for F is a quadrilateral having an 
elliptic vertex at z = i. In this case N-p{J-) has three elements. 

We associate to geodesic flow on Sj/T the (two-sided) cutting sequence subshift S-p,r which 
is a closed subshift of the full shift on |A'r('P)| letters. 

Definition 2.1. The set p of general position geodesies for V consists of those geodesies 7 
on ^ which intersect the interior of V, and which are transverse to all {gij^) : 5 G F} in the 
sense that 7 contain no vertex nor part of any side of positive length of any translated polygon 

^For non-convex V this may fail to hold. 
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Figure 2.1: The cutting-sequence encoding of a geodesic. The pictured geodesic has cutting 
sequence . . . , qq, gi, §2, ■ ■ starting with gi from V. 

{gV : g £ T}. If "P contains ideal vertices ("cusps"), then the endpoints of 7 must not coincide 
with any F-translate of any such vertex. 

To each 7 E G^^^, we assign the doubly-infinite cutting sequence 

C(7) := (• • ■ ,9-2,9-1,90,91,92, ■ ■ •) (2.2) 

in which ah gi G N^iV), as fohows. We first label each "side" s of the fundamental domain V 
with the label g{s, s') on its "inside" edge in V. We similarly label all the sides of the translated 
domain {hV : h £ T}, so that the edge of hV which is hs has the label g{s, s'). Every "side" is 
now assigned two labels, on its "inside" edge and "outside" edge, because it abuts two copies 
of a fundamental domain. (Figure 2.1 below indicates such labels.) The geodesic 7 = (^1,^2) 
visits in order a sequence of translates {hiV; z S Z} in which /ij G F and ho is the identity. At 
a crossing of the edge from hi-iV to hiV we assign the symbol gt on the side of the domain 
that the geodesic enters. (This is the convention of Katok [^^; and opposite to that of Adler 
and Flatto [3, p. 243], who use the symbol gi on the exit edge.) This convention yields: 

Lemma 2.1. The cutting sequence {gi, . . . , gj) follows a geodesic from V to the translated 
domain Vj = hjV with 

hj = gm- ■ ■ 9j- (2.3) 

Proof. The edge hiS of hiV that borders /ij-iP is assigned the edge label inside hiV of 
gi = g{s,s'). Now we have {hi^ig{s., s')h'^\)hi^iP = hiV, hence 

hi-igi = hi , (2.4) 

see Figure 2.1. I 

The transversality hypothesis says that if 7 intersects some translate g{V) for g GT, then 7 
intersects the interior of giV), which implies that the expansion ( |2.2| ) is well-defined and unique. 
If 7~^ = {O2, 9i) denotes the reversed geodesic, then G p has the cutting sequence 

C(7-^) = (... ,5-1, 50,51, •••) (2.5) 
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in which g'^ := {g_i) ^ 

Let p denote the set of all cutting sequences C(7) for 7 G Gpp. This set is invariant 
under the (forward) shift operator cr{{gi}) = {gi+i} since 

a(C7(7)) = C(5o7) , 

and the geodesic goj is in p. 

Definition 2.2. The orbit or shift- equivalence class [C] of a cutting sequence C is the union 
of all its forward and backward shifts 

[C] := U a'iCy, (2.6) 

i.e., it is the smallest shift-invariant set containing C. 

Thus [C(7)] is contained in S^, but is generally not a closed set. The geodesies on the 
surface S)/T are projections of geodesies on S). Since the shift operator on cutting sequences 
corresponds to a motion of the geodesic by an element of T, the orbit [C] is an invariant of the 
projected geodesic. 

Definition 2.3. The cutting sequence shift S-p^r is the closure of p in the symbol topology. 
That is, S-p^r consists of all symbol sequences 

,9-1,90,91,- ■ ■) 

such that every finite block {gi, gi^i, . . . , gi+k) occurs in some C{'~f) for a general position 
geodesic 7 € Gp p . 

We now extend the definition of cutting sequences to apply to all geodesies 7 that hit the 
interior of V. The set of cutting sequences C(7) for a general geodesic 7 = {6' , 9) consists of all 
symbol sequences that can be obtained as a limit point (in the sequence topology) of a sequence 
of G(7j) having 7j G p such that 7^ = (6*^, 6j) with 6'^ — > 6' and Oj^O as j^oo. There 
always exist such convergent sequences {jj : j > 1}, because Sp,r is compact in the symbol 
topology; hence G(7) 7^ 0. The set G(7) may be infinite for some geodesies 7. All cutting 
sequences in G(7) are contained in S-p^r; however, the closure operation used in defining S-p^r 
allows the possibility that S-p^r contains some symbol sequences not coming from any geodesic. 

Definition 2.4. The shift-equivalence class [G(7)] is the union of all orbits of cutting sequences 
in G(7). 

We show in the specific case F = PSL{2, Z) and fundamental domain that each [G(7)] is 



a union of a finite number of orbits of individual cutting sequences (Theorem 7.2). 

The totality of all possible convex polygons V that are fundamental domains of some group 
r has an explicit characterization. A general treatment appears in Maskit [^5|, which covers 
non-convex fundamental domains, and also covers groups of isometrics of Sj, which may include 
orientation-reversing isometrics. Maskit section 2] gives a sufficient condition for V to be 
a fundamental domain, but for finitely-generated groups where Sj/T has finite volume and V is 
hyperbolically convex it is a necessary condition as well. 
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.1 + v^i 
2^2'' 



S2 



g{si,s2) = [oi] = R 



g{s2,si) = [o 1 ] = L 



5(53,54) 



5(54,53) 



.2^2' 



1 
-1 



1 
-1 




Figure 2.2: The fundamental domain with side-pairings. 

2.2. Cutting Sequences for the Modular Surface 

We now speciaHze to cutting sequence expansions for geodesies on the modular surface 
Sj/PSL{2,'Z). The modular group PSL{2,7j) = SL{2,'Z)/ it I has many different convex fun- 
damental domains V that are hyperbolic quadrilaterals; see Beardon ||^, Example 9.4.4]. The 
standard fundamental domain for the modular group acting on f) is the hyperbolic triangle 



= {z : \z\ > 1 and 



-- < Re{z) < -} 
2 - ' - 2^ 



(2.7) 



which has one ideal vertex ("cusp") at ioo. When regarded as a fundamental domain, it is a 
quadrilateral having an elliptic vertex added at z = i, see Figure 2.2. 

The four "sides" of T are indicated in Figure 2.2 along with the two side pairings of s\ with 
S2 and S3 with 54, and with elements g{s,s') G PSL{2,7j), which comprise 



1 -1 
1 



1 1 
1 



1 
-1 
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The "sides" S3 and S4 containing the elhptic vertex together make up one side of J- ^ and the 
three sides are labeled inside the fundamental domain as follows: 



R = 



1 1 
1 



1 V3 



1 




1 



1 
-1 



S2 ■-- 



S3 U S4 : 




These matrices represent the PS'L(2, Z)-motions needed to move a neighboring fundamental 
domain J^' in to J^. These neighboring fundamental domains are labeled with the appropriate 
symbol in Figure 2.2. 

The action of PSL{2, Z) tiles S) with copies of J^, and the boundaries of all tiles fit together 
to cover a countable collection of geodesies of Sj, which are exactly those geodesies having two 
rational endpoints (f , |) such that 



det 



a c 
b d 



±2 , 



(2.9) 



where we adopt the convention that the cusp ioo is the rational ^. The resulting labeling of 
the tiles with the labels {L, R, J} is indicated in Figure 2.3. 

The cutting sequence shift := '^t^psl{2,'E) is defined by the method of section 2.1, and 
is a closed subshift of the shift on three letters {L,R,J}. For this special case, we present 
some additional cutting sequence labelings that keep track of corners, which we call generalized 
cutting sequences. We label the two finite corners of T symbolically by 



Ci = JRJ = LJL 
C2 = JLJ = RJR 



-1 

1 -1 

-1 

-1 -1 



1 V3 



(2.10) 
(2.11) 



General geodesies on the modular surface S^/ PSLi2,'L) are given by geodesies on ^ that 
pass through the fundamental domain. A representative set of lifts covering every geodesic in 
the modular surface consists of all vertical geodesies with —\<0<^, plus all semicircular 
geodesies in f) whose maximum point lies in the interior of the fundamental domain. The latter 
must necessarily hit the vertical line {it •.t> 0} and have endpoints (^1,^2) satisfying 

\02 - ^i| > V3, and - 1 < ^1 + 02 < 1. 

These conditions imply that 9i92 < 0. 

A generalized cutting sequence for a general geodesic 7 is the sequence of symbols from 
{L, R,J,Ci,C2} that describe the successive sides of the images of the fundamental domains 
that it hits. 



9 




Figure 2.3: Labeled tessellations of PSL{2,Z). The pictured geodesic, starting from x, has 
cutting sequence beginning R, R, J, L, L. 
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Finally, for the domain T we define a set of vertical cutting sequences which are one-sided 
cutting sequences associated to the set of oriented geodesies of that emanate from the cusp 
{icxD} of T and pass through T , which are exactly the set of vertical geodesies {(cx),^) : — ^ < 



of those vertical geodesies {oo,9) for irrational 6 with — i < < ^ 



Definition 2.5. The irrational vertical cutting sequence set 11^ consists of cutting sequences 

l<0<i. 



An important result for our analysis is the following simple fact. 

Lemma 2.2. An irrational vertical geodesic cannot hit any finite corner of a PSL{2,7j)- 
translate of J- . 

Proof. The generators of PSL{2, Z) acting on the upper half-plane are z ^ z + 1 and 
z — > —1/z. Both generators preserve the field Q(\/— 3), in which every element has rational 
real part. The finite corners ±1/2 + -v/— 3/2, and thus all of their images, are in this field. I 

Thus the cutting sequences in 11^ are one-sided infinite sequences which contain no symbol 
Ci or C2, hence they are all of the form: 

C+W := (50,51,92, •••) all 5iGiV(.F) = {L,R,J} . (2.12) 

The expansion C~^{9) is a kind of continued fraction of 0, and we call it the cutting sequence 
expansion of 9. Note that, for a given point 6 + it corresponding to a fundamental domain hiJ-, 
that there is some point z ^ T such that 

9o9i ■ ■ ■ 9i-29i-iiz) = 9 + it, (2.13) 



according to Lemma 2.1 



Definition 2.6. The vertical cutting sequence space Iljr is the closure of in the one-sided 
sequence topology on the one-sided shift space on three symbols {L, R, J}. 

Note that the set 11^ is not closed under the one-sided shift operator o""*", because the first 
symbol of every element of 11^ is J. 

We can now define cutting sequences C{9) for rational 9 with —■^<9< ^ in IIjp as limits 
of cutting sequences C{9i) of irrational 9i — > 9, as described earlier; in this case C{9) is a set of 
one-sided infinite cutting sequences. 

We define generalized cutting sequence expansions for vertical geodesies, drawn from the 
alphabet {L, R, J,Ci,C2}, similarly to the case of general geodesies. 



3. Continued Fraction Expansions and Cutting Sequence Expansions 

This section describes symbolic dynamics for vertical cutting sequences associated to ad- 
ditive continued fraction (ACF) expansions, and for Minkowski geodesic continued fraction 
expansions (MGCF) for vertical geodesies. It shows that the cutting sequence expansion and 
MGCF expansions for a geodesic can each be computed from the other using a finite automaton. 
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3.1. Additive Continued Fraction Expansions 



The ordinary continued fraction (OCF) expansion for a real 9, written 9 = [ao,ai,a2, ■ ■ •], 
can be represented using 2x2 matrices. The use of such matrices to represent continued 
fractions appears in Frame |13] and Kolden |28], and is described in Stark [51|. Namely, one 
has 

ao 1 1 \ an ^ \ Pn Pn-l 1 /o iN 

.1 J ■ • • [ 1 J [qn Qn-l J ' ^ ' ' 

where |^ = [ao,ai, • • • ,an] is the n-th convergent of 9. The ordinary continued fraction expan- 
sion for 9 > has a symbolic dynamics which is a one-sided shift on infinitely many symbols, 
using the alphabet 



k 1 
1 



A; > 0. 



(3.2) 



The additive continued fraction (ACF) expansion of > is a symbolic expansion using the 
two symbols 

"ill r n 1 " 

(3.3) 



(3.4) 





1 ■ 




" 


1 " 






, F = 







1 




1 






which is obtained from (|3.2[) by expanding the symbol Lj^ as 



' k 1' 




" 1 


1 ■ 


k 


"01" 


10 







1 




10 



The additive continued fraction expansion of ^ > is thus 

ppaippaappaa ... 



(3.5) 



regarded as a string of letters in the alphabet {R, F}. The finite truncations of this expression 
are 2x2 matrices that encode the convergents and intermediate convergents of 9. 

We next construct a related continued fraction for real 9 > which uses the symbols 



R 



1 1 
1 



D 



1 
1 1 



(3.6) 



The Farey tree expansion of a real number > that has OCF expansion 9 

paoQaipa2Qa3 



[ao,ai,a2, 



is: 



(3.7) 



regarded as a sequence of letters in the alphabet {R,D}. This sequence encodes the steps of 
subtraction needed in the division process to encode the ordinary continued fraction of 9] it 
also essentially gives all the intermediate convergents to 9, cf. Richards [Q], Theorem 2.1. For 
^ > the set of allowable symbol sequences for the additive expansion is the full one-sided shift 
on two letters {R, D}. 

A relation of this expansion to paths in the Farey tree is described in Lagarias [^] and 



Lagarias and Tresser [32|. We start with the interval [y, represented as the matrix 



1 

1 



At each step, we multiply our current matrix 



P2 q2 



on the left by R or D, as appropriate. 
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The matrix R causes our current matrix 
matrix D causes it to be replaced by 



to be replaced by 



Pi +P2 qi + q2 

P2 q2 



; the 



, Thus, in either case, the next term 



Pi Qi 
P2 q2 
Pi qi 
Pi +P2 qi + q2 

in the Farey tree is {pi +P2)/{qi + Q'2)) and we continue with the interval on the side of this 
approximation which contains 9. 

The symbolic expansions ( |3.5D and (^]^) are equivalent in the sense that each is convertible 
to the other using a finite automaton; see section | 

3.2. Cutting Sequences and Minkowski lattice basis reduction 

We now develop a correspondence between the cutting sequence expansion of the geodesic 9 + it 
for — ^ < 9 < ^ and the construction of a Minkowski-reduced lattice basis for the parametrized 
series of lattice bases 

1 



t 



(3.8) 



for the parametrized family of lattices At = 0), {—9, t)]. 

To demonstrate the equivalence, we first transform lattice bases in GL(2, M) to the cone of 
positive definite symmetric matrices V2 (equivalently, to positive definite quadratic forms), and 
then we relate these to elements of Sj. 

Definition 3.1. A lattice basis vi,V2 is Minkowski-reduced if vi is a shortest vector in the 
lattice, and V2 is a shortest vector not a multiple of vi. 

For each t, there is a matrix P = (possibly several) for which PB((^) is Minkowski- 
reduced. General properties of Minkowski-reduction can be found in Cassels [|lO| and Gruber 
and Lekkerkerker |]l5|, pp. 149ff]. 

A basis B G GL(2,M) corresponds to the positive definite symmetric matrix M = BB"^ G 1^2 
with associated quadratic form 



Xl 
X2 



mnx? + (mi2 + m2i)xiX2 + m22xi 



Q{xi,X2) = [xiX2]M 

For the basis matrix Bf(0) in ( |3.8D this gives 

M := Mt{9) = 

with associated quadratic form 

Qt{xi,X2) = {xi - 9x2)"^ + t^x\ = (xi - {9 + it)yi){xi - {9 - it)x2) 



(3.9) 



(3.10) 



(3.11) 



The action of P G GL(2,Z) which sends the matrix ^t{9) to PBt(0) is transformed to the 
GL(2, Z)-action on 1^2 that takes M to PMP"^, and this takes the associated quadratic form 
Q{xi,X2) to 

Q'(xi,X2) = Q{piixi +P21X2, P21X1 +P22X2) ■ (3.12) 



13 



The Minkowski reduction domain for GL{2,Z,)\V2 (see Cassels pO[] ) is the set of quadratic 
forms (3^) such that 



g(i,o) < Q(o,i), 

Q(0,1) < Q(l,l), 
Q(0,i) < Q(-i,i). 

This is equivalent to the conditions 

\mi2 + m2i\ < mil < 



(3.13) 



(3.14) 



on the form Q in (O). The subgroup 

-'-{111 



-1 
1 



1 

-1 



-1 
-1 



of GL(2,Z) maps the domain M into itself, hence we actually let H2\GL(2,Z) act on V2- 
By definition the Minkowski reduction domain Ai of GL(2, Z)\GL(2, M), is the preimage of 
A4 under the mapping B BB^, and the group actions by GL{2, Z) coincide. An algorithm 
for choosing Minkowski-reduced lattice bases will choose one element of GL[2,'L) out of the 
four-element coset of H2\GL{2,Z). A natural canonical choice is to require det P = 1; this 
still leaves two choices for P, but multiplying P by — 1 will not change the quadratic form and 
can be done arbitrarily. A different choice will be made by the Minkowski geodesic continued 
fraction of the next section. 

For the second transformation from 7^ to ^, note that a positive definite quadratic form can 
be uniquely written 

Q{x,y) = a{xi - 9x2){xi - 9x2) , (3.15) 

with 9 and 9 complex conjugates and Im{9) > 0. We map P to ^ by sending the quadratic form 
Q to the unordered pair of roots {9, 9} and then identify this with 9 £ S) = {z : Im{z) > 0}. If 
P G GL(2,Z) sends Q to Q' given by ( |3.14| ), then a calculation yields 

Q' = a'{xi - <pX2){xi - 4>X2) , 



with a' = Q{pii,Pi2) and with a root (j) given by 



P22 
-Pi2 



-P21 
Pu 



(3.16) 



where the action of the matrix is a linear fractional transformation on C. If det(P) = 1, then 
4> € Sj, while if det(P) = —1 then Im{(p) < is the complex conjugate of the root we want, i.e. 
it reverses the ordering of the roots. Now each four element coset of H2\GL{2, Z) contains two 
elements in SL{2,'L) which form a coset of PS'L(2,Z) = {ibl}\S'L(2, Z), so the action ( 3.16| ) 
of i?2\G'L(2,Z) on unordered pairs {9,9} is equivalent to the standard PS'L(2, Z)-action on i^. 
Finally, the conditions ( 3.13| ) for a form Q to be Minkowski reduced translate via ( |3.15| ) to 



> 1, 



-1>9 + 9>1 



which is exactly the fundamental domain T of PSL{2, 
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The mapping from GL(2,M) to Sj obtained by composing these two transformations sends 
the set of matrices {Bt{G) : t > 0} in GL{2,R) to the geodesic {d + it : t > 0} in S). Since the 
H2\GL{2,'L) action on GL(2,M) is equivalent to the PS'L(2, Z) action on ^ and the reduction 
domains correspond, the sequence of Minkowski-reduced lattice bases for 6 is essentially the 
same as the cutting sequence expansion for 0. In the reverse direction, 6 €z Sj determines a 
positive ray {aQ : a > 0} in the cone V2, i-e. an element of 1^2/^^^, and this in term determines 
a family {aBQ : a G R+,Q G 0(n,M)} in GL(2,R), i.e. an element of GL(2,R)/R+0(2,M). 
This does not affect the symbolic dynamics because the Minkowski domain M is invariant 
under the M+-action, and because the Minkowski reduction domain M is invariant under the 
i?*0(2,M)-action. 



3.3. Minkowski Geodesic Continued Fraction Expansions 

The Minkowski geodesic multidimensional continued fraction (MGCF) expansion is intro- 
duced in Lagarias [30|. We consider here the one-dimensional case. 

This is a specific algorithm for Minkowski-reduction of the lattice bases 



1 
-9 t 



(3.17) 



There exists a sequence of critical values 



00 



to> ti> t2> 



with lim„^oo tn = 0, and an associated sequence of convergent matrices {P^") : n = 0, 1, 2, • • •} 
in GL(2,Z) such that P'-"'^Bj(0) is a Minkowski-reduced lattice basis of the lattice when t„ > 
t > tn+i- We write 

r {") 



(n) 

(n) 
Q2 



(3.18) 



so that 



p(")Bt(0) 



(n) 
Pi 

(n) 
P2 



t 



(3.19) 



This motivates the name "continued fraction" ; as t goes to zero, the approximations p\ 
must also go to zero, and thus the pi/qi play the role of convergents. 

Changing the sign of one or both basis vectors does not affect Minkowski-reduction, so we 
have four choices of P*-"^. For defining the MGCF as a continued fraction, it is natural to 

(n) (n) 

require positive denominators; that is, we require that q"^ ' > for ah n > I, and if q); ' = 
then p-"^ > 0. For lattice basis reduction, it is natural to require det P^*^-* = 1; multiplying by 
±1 is still possible but will not affect the process. 

In either case, the associated partial quotient matrices A^"^ are defined by 

p(n) _ ^(n)p(n-l) 



SO that 



p(n) ^ p^in) i^{n-l) . . . ^(l)p(O) 



(3.20) 
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Here P^'^^ = I is the identity matrix if — ^ < ^ < ^, since the lattice basis {{1,0), {—9, t)} is 
Minkowski-reduced for large t. 

We now show that the partial quotient matrices A^"^ are drawn from a finite set. 

Lemma 3.1 (Minkowski partial quotient set) . The allowed partial quotients for the one- 
dimensional Minkowski geodesic continued fraction algorithm are 



1 

1 



1 
1 1 







1 1 
1 



(3.21) 



Proof. Let M 



vi 

V2 



be a 2 X 2 matrix with det(M) > 0, and associate to it the positive 



definite quadratic form with coefficient matrix 



|vi| 



(vi, V2) 



(vi,V2) ||v2r 

The Minkowski-reduction conditions for the quadratic form Q{x.) = x^(MM^)x are 



|vi| 

|V2| 
|V2| 



< 
< 
< 



|V2|| 

|vi + V2I 
|vi - V2I 



(3.22) 



(3.23) 
(3.24) 
(3.25) 



see Cassels |1C, p. 257]. 

We must now choose our algorithm; we will describe both the natural algorithm for Min- 
kowski lattice basis reduction and the Minkowski geodesic continued fraction algorithm. 

For lattice basis reduction, it is natural to apply a new matrix A*-"^ of determinant 1. At 
any critical time t„, all of the above inequalities are satisfied with the current matrix p("^^) 
for tn-i > t > tn, and at least one holds with equality at t — t^, 

We first consider "generic" convergents which occur when exactly two of ||vi(t)||, ||v2(t)||, 
||vi(i) -|-V2(t)|| and ||vi(t) — V2(t)|| become equal at t = tn- 

If at t = tn we have ||vi(t)|| = ||v2(t)|| then the partial quotient matrix is J = '^x 

it exchanges vi and V2 and changes the sign of the new V2 to keep determinant 1. If ||v2 

|vi(t) + V2(t)|| then the partial quotient matrix is R 1 1 



1 1 



If ||V2 



|vi(t) — V2(t)|| then the partial quotient matrix is L 



; it replaces V2 with vi + V2. 

; it replaces V2 



1 
-1 1 



with V2 — vi. 

We might also have more than one of the inequalities ( p.23D -( ^.25 ) holding with equality at 
the same t. If ||vi(t) + V2(t)|| = ||vi(t) — V2(t)||, then neither one can be equal to ||vi(t)|| unless 
V2(i) = 0, nor to ||v2(t)|| unless vi(t) = 0. We could have ||vi(t)|| = ||v2(t)|| = ||vi(t) + V2(i)||. 
In that case, for t < tn, we would have ||vi(f)|| > ||v2(t)|| > ||vi(t) + V2(t)||, and to correct 
these inequalities, we must replace vi by vi + V2; we can take partial quotient matrix Ci = 
1 —1 1 

by also changing signs of both vectors Similarly, if ||vi(t)|| = ||v2(i)|| = 



JLJ = 

l|vi(t) 



-1 



V2(t)||, we get partial quotient matrix C2 = JRJ 



-1 1 
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These five matrices (generated from just three) give the natural algorithm for Minkowski 
lattice basis reduction. 

For the Minkowski geodesic continued fraction, we follow a similar process except that we 
choose signs to keep the qi positive. We use the fact that all vectors in the lattice Aj(0) have 
the special form 

Vi = Vj(t) := (pi - QiO, Qit) with > 0. (3.26) 
Whenever a critical value t„ occurs with || (p — qO, qtn)\\ = \\ip' — q'S, q'tc)\\, then the inequality 



l't)\\ > 



,qt)\\ for 0<t<tc 



holds exactly when \q'\ > \q\. Thus when the shortest vector in the basis P^""-*^) is replaced by 
a vector in P^"^ the associated denominator must increase. 
Again, we first consider "generic" convergents. 

Suppose first that qi < q2- If at t = t„ we have ||vi(t)|| = ||v2(t)|| then the partial quotient 
; it exchanges vi and V2. If ||v2(t)|| = ||vi(t) + V2(t)|| then the partial quotient is 



|vi(t) — V2(t)|| cannot occur. 



it replaces V2 with vi + V2. The case ||v2 



1 
1 

1 1 

Suppose next that > (72- At t = t„ we cannot have ||vi(t)|| = ||v2(t)|| since an exchange of 
vi and V2 would not increase the denominator. We may have ||vi(t) + V2(t)|| = ||v2(t)||, giving 

1 



the partial quotient 




1 1 



If qi > 2q2 then q' = qi — q2 > <?2 hence the partial quotient 



is possible; it replaces V2 with vi — V2. 



1 

1 -1 

There remain situations in which three or more of ||vi(t)||, ||v2(t)||, ||vi(t) + V2(t)|| and 
vi(t) — V2(t)|| simultaneously becoming equal at t = tc- Only one case is possible; it is 



= I|v2(t)|| = ||vi(t)+V2(t)|| , 

which can only occur when qi < q2. In this case the partial quotient is 



The other 



1 1 
1 

case consistent with an increasing denominator is ||v2(t)|| = ||vi(t) + V2(t)|| = ||vi(t) — V2( 
but this implies ||vi(t)|| = 0, which is impossible unless t = 0. M 

We define the Minkowski geodesic continued fraction expansion of any real satisfying 
— ^ < ^ < ^ to be the symbol sequence 



(A(i),a(2),A(3),...) . 



(3.27) 



We write this sequence left-to-right, although the matrix product in ( p. 201 ) runs right-to-left. 
The associated Minkowski geodesic symbol set is 



L 



" 1 ■ 


,R = 


"10' 




" 1 " 




" 1 1 " 


,J = 


,c = 


1 


1 -1 




1 1 




10 





(3.28) 



The symbol C occurs only for "exceptional" geodesies, and we will usually be concerned with 
symbolic expansions drawn from the smaller symbol set {L, R, J}. 
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The Minkowski geodesic continued fraction expansion corresponding to a general geodesic 
on T is the expansion attached to Minkowski reduction of the parametrized lattice bases 



The associated quadratic form is 



1 {QxY^t 

-02 t 



1 



Qt{x, y) = {x- d2yY + t'[—x-y 



(3.29) 



and the partial quotients are obtained by the same formula. 

3.4. Correspondence between the continued fraction and cutting sequence 

We will now give a correspondence between the symbol sequences for given by cutting se- 
quences, the natural Minkowski lattice basis reduction algorithm, and the Minkowski geodesic 
continued fraction algorithm. 

The precise correspondence between the symbol sequences involves specifying the relation 
among the four elements of the i?2-coset in H2\GL{2,Z) for the Minkowski geodesic continued 
fraction expansion and the two elements in the {ibl}-coset in PSL{2, Z) for the cutting sequence 
expansion. 

Definition 3.2. The parity of a word W = Si . . . S^, in the alphabet {L, R, J} is even or odd 
according to whether there are an even or odd number of L and J. That is. 



det(W) 



-1 



,parity(W) 



(3.30) 



We obtain: 



Theorem 3.1. For irrational 9 with — ^ < 9 < ^, the one-sided cutting sequence expansion of 
the vertical geodesic (cxd, 9) = {9 + it : t > 0} in 11^ is obtained from the Minkowski geodesic 
continued fraction expansion of 9 by the following procedure: If the current initial word of the 
MGCF expansion has even parity, on the next symbol make the letter replacement L — > L, R — > R 
and J ^ J; if it has odd parity, make the letter replacement L — > R, R — > L and J — > J. The 
MGCF can be obtained from the cutting sequence by the reverse process, in which the parity is 
determined by the current symbols of the MGGF expansion. 



Proof. By Lemma 2^, the assumption of irrational 9 ensures that the MGCF expansions 
and cutting sequence expansions of 9 are both infinite and never use a symbol C. The MGCF 
expansion • • • uses the symbols 



L 



drawn from GL{2, Z) and moves right to left, while the cutting sequence expansion uses the 
symbols 



" 1 ■ 


,R = 


" 1 " 




" 1 " 


, J = 


1 -1 


1 1 




1 



1 -1 

1 



1 1 

1 



1 
-1 
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drawn from S'L(2,Z) and moves left to right by Lemma 2.1 
introduce the matrix 

' 1 



For notational convenience, we 



K 







-1 



so H2 is {l,K, — K, — I}. We first arrange for matrices to multiply in the same direction. The 
matrices in the MGCF are chosen so that P'^"''^ Bt{9) G ^ for a given t, while the matrices 
in the cutting sequence are chosen so that the cutting sequence product /i„ = Si . . . S„ has 
7(t) G hnJ-. Therefore, we have /i„ = (P^"^)~^. We expand (p("))~i in terms of symbols 
L, R, J, K using the relations = L, = KRK and = J to obtain an expansion which 
proceeds left to right. To convert this expansion to the cutting sequence expansion, we must 
convert to symbols L, R, J and remove the symbols K, which encode the parity. We first convert 



to the symbols of the natural Minkowski basis reduction algorithm of Section 3.3 



" 1 ■ 




" 1 " 


-1 1 




-1 



R := R, 



in SL{2,Z). To convert a product of matrices from the form L, R,J to L, R,J, starting from 
the right end of the product, we use the relations LK = R, RK = KL, and JK = —J. In doing 
this we pick up or lose a factor of K whenever we encounter a matrix L or J and this multiplies 
the determinant by —1; meanwhile, every R in the MGCF becomes R~^ = KRK when inverted, 
and thus is itself encoded with the opposite parity but leaves the parity unchanged for the next 
matrix. (We sometimes pick up a matrix factor of —I, but this commutes with everything and 
may be ignored.) Similarly, to reverse this, we use L = KL and J = KJ. 

To convert from the natural Minkowski basis reduction algorithm symbols L, R,J to the 

_ _ _ ~X' — " T — "T — 

cutting sequence symbols L, R, J, we use the relations (L = R, (R = L and (J = J. 
Thus the conversion from L, R, J to L, R, J interchanges R with L. I 



3.5. Finite Automata 

By a finite automaton we mean a deterministic finite-state automaton, as defined in Hopcroft 
and Ullman |^ or |43], used as a transducer. Such an automaton is a finite directed graph 



with labeled edges, which may contain loops and several edges exiting from each vertex. The 
states are the vertices of the graph, and the edges give rules to move from one state to the 
next. Each edge has two labels, an input label and an output label. If the symbol alphabet has 
s letters, then from each vertex there exit exactly s edges whose input labels are exactly the s 
allowed symbols. The output labels are finite strings of symbols, possibly empty. The machine 
starts in a given state. It reads an input symbol, which tells it which exit edge to follow, prints 
the specified output string of letters as output, and moves to the new state specified by the 
edge. Then it proceeds to the next input symbol. 

In this paper we present a number of results asserting the existence of finite automata to 
convert one-sided infinite symbol sequences of one form to another form. In the proofs we only 
indicate the "finite-state" character of the conversion process, and generally omit details of the 
(routine but sometimes involved) construction of the automaton. 

As a simple example, the discussion in section 3.1 yields: 

Theorem 3.2. For real 9 > 1, the additive continued fraction expansion of 9 can be converted 
to the Farey tree expansion of 9 by a finite automaton, and vice versa. 
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F:0 



(a) Additive continued fraction to Farey tree 



R : R 



J : J 



R : L 




(b) MGCF to cutting sequence 
Figure 3.1: Finite Automata-Transducers (initial state for both automata is +1). 



Proof. To convert from the Farey shift expansion (^J) to the additive continued fraction 
expansion (^.Sj), we use D = FRF and F^ = I. 

A finite automaton which converts the additive continued fraction expansion ( |3.5D to the 
Farey tree expansion ( |3.7D must keep two states to keep track of the sign ±1 of det(Si • • • Sm) 
of the symbols Sj = F or R examined so far. The initial state is +1. I 

The discussion in section 3.4 yields: 

Theorem 3.3. For irrational 9 with — | < 9 < ^ there exists a finite automaton to convert 
the Minkowski geodesic continued fraction expansion of 9 to that of the vertical cutting sequence 
expansion of (oo, 9) and vice-versa. 



Proof. This follows from Theorem 3.1 



For each direction, the finite automaton constructed 
needs two states, to keep track of whether there a factor of K present; i.e., to keep track of the 
sign of det(Si • • • Sn). The initial state is +1. I 

We illustrate two of these automata in Figure 3.1. The edge symbol S : W specifies that 
this edge is taken if the current input symbol is S, and W denotes a symbol sequence to be 
output. The initial states are the vertices labeled +1. 



4. Vertical Cutting Sequences to Additive Continued Fractions 

The ordinary continued fraction expansion of an irrational real number — ^ < ^ < ^ is easily 
determined from either its cutting sequence expansion or its Minkowski geodesic continued 
fraction expansion. We begin with the latter case. 
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Theorem 4.1. Given 6 with — 2 < ^ < 2? with ordinary continued fraction expansion 9 = 
[oo; ai, a2, 03, • • •], let S0S1S2 • • • be the Minkowski geodesic continued fraction expansion of 9 in 
the alphabet {L, R, J, C}. The MGCF expansion can be uniquely factored into segments BqBi • • • 
where Bq = J or JL and each Bj for i > 1 is R^J for some k > 1, R'^'''^JL for some k > 1, or 
R^C for some k > 1. Each segment encodes exactly one or two partial quotients of the OCF 
expansion. For general segments, R^J encodes a partial quotient an = k, while R'^^^'^JL or R^C 
encode two partial quotients an = k, a„+i = 1. For the segment Bq, the symbol J encodes 
oq = 0, while JL encodes the partial quotients ao = —1, ai = 1- The case R^C can only occur 
when 9 is rational. 



Proof. We say that p/q is a best approximation to 9 i( \q9 — p\ < \q'9 — p'\ for < < q. 
Using this definition, it is a basic result about ordinary continued fractions that the convergents 



are the complete set of best approximations; see Hardy and Wright |17, Theorem 182]. We say 
that p / q is a better approximation to 9 than p' / q' \q9 — p\ < \q'9—p'\. The best approximation 
property thus proves that Pn/Qn is a better approximation than any fraction with denominator 
less than qn+i- 



We start with P^^^ 



the identity matrix. Thus P^^^Bt{9) is 



1 





t 



Neither 



of the Minkowski inequalities ( 3.2^ ) and ( 3.2E| ) can hold with equality for any t, since \9\ < ^. 



We get ti = 1 — [\_9 + ^\ — 9) as the value at which ( p.23|) holds with equality. The first partial 



quotient matrix is thus J. That makes P^^^ 




1 



. Now P^^^Bt{9) 



1 



t 




We check the Minkowski inequalities. Inequality ( p. 231 ) is satisfied, with equality only at 
t = ti. Inequality ( 3.24 ) is an equality ifl = (1 — ^)^+t2) and ( ^.251) is an equality if 



1 = (—1 — 9) + ^2- If ^ = 0, neither one holds with equality for any t > 0, and J is thus the 



whole MGCF, encoding 9 = [0]. If < 6* < i, then only (|3^ ) can hold with equality, so the 
next partial quotient matrix will be R. In this case, we have 



p(i) 



" 


1 ■ 




PQ 




1 









q-i _ 



and we let Bo = J. If -i < 6* < 0, then only ( pSj ) can hold with equality, so the next partial 
quotient matrix will be L. For —^<9<0, the continued fraction for 9 begins [—1, 1, . . . ], and 
thus pi = 0,qi = 1. This gives 



p(2) 



and we let Bq = J, L. 

In either case, after the first step, we have encoded all coefficients through a„_i, and our 
current matrix is 



" 1 ■ 




Pi qi 


-1 1 




_ po qo _ 



p(i) 



Pn-l 
Pn-2 



ln-1 
ln-2 



(4.1) 



The rest of the theorem is proved by induction, with the induction hypothesis that after each 
segment B^,, the matrix P^-'^ is in the form ([4.1|), with the coefficients through a„_i encoded. 
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First, assume that the current P^-^^ is in the form (4.1). Then 



Pn-l -Qn-lO tqn-l 
Pn-2 — 9n-2^ tqn-2 



Since Qn-i > Qn-2, the requirement of an increasing denominator makes J impossible for the 
next partial quotient. The Minkowski inequality for R is 



{{Pn-l - Qn-lO) + {Pn-2 " qn-20)) - {Pn-2 " 



+ t^{{qn-i + qn-2f - q, 



n-2) 



> 0, 



(4.2) 



and for L, it is 



Hp 



'n-l 



Qn-lO) - {Pn-2 " qn-20)f - {Pn~2 " qn-20f 



dn~2) 



> 0. 



(4.3) 



By the best approximation condition, pn-i/qn-i is a better approximation to 6 than (pn-i — 
Pn~2)/ (qn-i — qn-2)- Thus ( [4. 3D is Satisfied for small t, and thus for all t < tj, so we cannot have 
L next. If 9 is rational and fln-i is the last term, then ( [4. 21) is satisfied for all t; thus the MGCF 
terminates. Otherwise, since every even convergent is less than 9, and every odd convergent is 
greater, Pn-i — qn-i9 and Pn-2 — qn-29 have opposite signs, and the absolute value of their sum 
is less than the absolute value of Pn-2 — qn-29- Thus ( [4. 2D is not satisfied for sufficiently small 
t, and tj+i is the value of t at which it holds with equality. Thus the next matrix is R; hence 



Pn-l qn-i 

Pn-2 + Pn-l qn-2 + qn-1 



p(i+i) _ Rp(i) 

The remainder of the proof considers the case where P^-'-' is of the form 
p(i) 

which can be reached after m applications of R to the form ( |4.1| ). Now we have 
P^'^Bt{9) 



Pn-l qn-1 
Pn-2 + mpn-l qn-2 + mqn-1 



(4.4) 



Pn-l -qn-l9 tqn-l 
{Pn-2 - qn-29) + m {pn-1 - qn-l9) t{q 

n—2 1 Tnqn—1) 



Since qn-2 + mqn-i > qn-i, the next symbol yl(-^+^) cannot be L. The other two Minkowski 
inequalities, for J and R, are 



2 2 
i{Pn-2 - qn-29) + m {pn-1 - qn-l9)) - {Pn-1 " qn-l9) 

+ t\{qn-2 + rnqn-i)"^ - ql^^) > 0, 

((m + 1) {pn-1 - qn-l9) + {Pn-2 " qn-29)f 
- {m {pn-l - qn-l9) + {Pn-2 " qn-29))^ 

+ t^{{{m+l)qn-i + qn-2)'^ - {mqn-1 + qn-2)'^) > 0. 



(4.5) 



(4.6) 
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Thus the next symbol is A'--'^^^ = J if (}4.5D holds with equality for a larger t < tj than ([4.6|), 
A(i+i) = R if holds with equality for a larger t < t^, and A(^'+^) = C if both hold with 
equality at the same t. Recall that we have, for < m < a„, and n even, 



Pn-l + Pn {an + IjPn-l + Pn-2 ^ ^ '^Pn-l + Pn-2 
~ > a > 



Jn-1 + Qn {an + l)qn-l + qn-2 



mqn-1 + qn-2 



(4.7) 



with the reverse inequalities holding for n odd. The last inequality holds with equality only if 
9 = Pn/qn and m = an- If m < a„, then mqn-i + qn-2 < qn, so we must have 

{pn-l - qn-lOf < {{Pn-2 " qn-20) + m {pn-l - qn-lO)f, 

and thus the left-hand side of (|4.5| ) is positive for small t, so A'--'^^-* cannot be J or C. Thus, 
if m < Qn, the only possible case is R, and the matrix p(-'+^) is still of the form ( [4. 41 ), with m 
replaced by m + 1 . 

We next consider the case when m = in ([4.4| ). In this case, the next symbol A^-'^"'^^ may 
be any of R, J, and C. First suppose A^-'"'"^^ = J. If so, we have finished our segment = R"" J, 
and P^-'"*"^) is now in the correct form (14 .Jj), since it is 

Pn qn 

Pn-l qn-1 

Next suppose that m = an and A^-''*'^^ = R. For this to happen, we need ( [4.6|) to be true 
only for t > t^. Since m = an, we have mqn-i + qn-2 = qn, and similarly for p„. For ( |4.6D to be 
true only for sufficiently large t, {pn-i + Pn) / {qn-i + qn) must be a better approximation to 6 
than Pn/qn- This cannot happen ifO = Pn/qn, hence a„_|_i is defined. By the best approximation 
property, we must have qn-i + qn ^ qn+i, and it follows that ffln+i = l, Pn+i = Pn + Pn-i, and 
qn+i = qn + qn-i, SO that in this case 



p(j+i) _ Rp{i) 



Pn-l qn-1 
Pn+1 qn+1 



This matrix is in the form ( [4.4| ) but with m = a„ + 1. Now A'--^"''^-* cannot be R or C, because 
the left-hand side of ( |4.6| ) is positive for small t, so must be J, which gives 



p(i+2) ^ JRpO') 



Pn+1 qn+1 
Pn-l qn-1 



(4.8) 



This gives 



p(^'+2)Bt(0) 



Pn+1 — qn+lt 
Pn-l - qn-li 



tqn+1 
tqn-1 



The following symbol A'--'"'''^^ cannot be J, because qn+i > qn-i and we must have an increasing 
denominator. The Minkowski inequality for R is 



{{Pn+l - qn+lO) + {Pn-1 " ^n-l^*)) " {Pn-1 " qn-lO) 

+ t^{{qn+i + qn-l)^ -ql-l) > 0. 



(4.9) 



23 



The constant term here is positive because Pn+i — Qn+iO and Pn-i — Qn-iO have the same sign, 
so this inequaUty holds for all t. The Minkowski inequality for L is 



{{Pn+l - Qn+lO) - (Pn-l " Qn-lO)) - {Pn-1 " Qn-lOY 



In-l) 



which is not satisfied for small enough t because {pn+i — Qn+iO) — {Pn-i 
and thus its constant term is negative. Thus A^-''^^^ = L, so that 



> 0, 



p(j+3) = LJRP(^') 



Pn+l Qn+1 
Pn Qn 



(4.10) 

Pn - QnO, 

(4.11) 



This matrix has the form ( [4.1[ ), and we have encoded the two coefficients, a„ and fln+i = 1, 
with a segment B^. = R""+ML. 

Finally, suppose that m = and A(^+i) = C. a gain, ( |4.6| ) must fail to hold for sufficiently 
small t, and it follows that a„+i = 1. Thus, since Qn+i = Qn + Qn-i and Pn+i = Pn + Pn-i, we 
have 



p(j + l) ^ Qp(j) 



1 1 

1 



Pn-1 Qn-l 
Pn Qn 



Pn+l Qn+1 
Pn Qn 



(4.12) 



Here p(-'+^) is in the form ([4.1|) , and we have again encoded the two coefficients, a„ and a^+i = 1 
as R""C. This completes the induction step. 

Lemma \i.2\ shows that the last case R'^C can occur only for rational 9. I 
There is an analogous conversion method from the cutting sequence expansion to the addi- 
tive ordinary continued fraction expansion, as follows. 

Theorem 4.2. Given with ordinary continued function expansion 6 = [oq; oi, Q^S; ■ " 
808^82 • • • be the cutting sequence expansion for the geodesic (00, 9) = {6 + it : t > 0} in the 
alphabet {L, R, J, Ci, C2}. It can be uniquely factored into segments 606162- •• where 60 = J 
or JR, encoding oq = or oq = —1, ai = 1, respectively, and each succeeding segment is 
R^J or L^J encoding an for n even and odd, respectively, or is r'^^'^J R, or R^Ci encoding 
a„ = fc, a„ + 1 = 1 for n even, or is l'^^^J L or L^C2 encoding an = k, a„+i = 1 for n odd. The 
symbols Ci, C2 can only occur in expansions of rational 9. 



Proof. This follows from Theorem [4.1| by noticing that the parity of the initial word 
Wj = 60 ... 6, changes after each segment 6j encoding one term, and does not change after any 
segment 6j encoding two terms, while the segment 69 has odd parity if it encodes no terms, 
and even parity if it encodes ai. I 



An immediate consequence of Theorem 4.1 and 4.2 is the following result 



Theorem 4.3. There exists a finite automaton which converts the Minkowski geodesic contin- 



ued fraction expansion of each irrational 9 with 



<9 < 



to the additive continued fraction 



expansion of 9. There exists a finite automaton that converts the vertical cutting sequence 
expansion of {00, 9) to the additive ordinary continued fraction expansion of 9. 



Proof. The segment-partition of the MGCF expansion given in Theorem |4.1| permits the 
Farey tree expansion (|3.7| ) to be computed by a finite automaton, because the necessary infor- 
mation to decide on the symbol R versus L depends only on the determinant ±1 of the product 
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of the MGCF matrices scanned plus the values of the last four MGCF symbols in the expansion. 



Next, Theorem [4.2| guarantees that a finite automaton exists to convert the cutting sequence 



expansion as well. Finally, Theorem 3.2 applies to give the additive continued fraction from 
the Farey tree expansion. I 

Theorem [4.1| implies that the Minkowski geodesic continued fraction expansion of 6 can 
represented in an abbreviated form 

6 := [00,01,02,03, ...] 
similar to its ordinary continued fraction expansion 

9 = [00; 01, 02, 03, ...] , 

with the change that each symbol 1 in the OCF expansion is to be replaced by one of three 
possible symbols 1^, Im and Ic- Here Ih means that the continued fraction partial quotient 
On = 1 begins a new segment RJ or RRJL in the MGCF expansion (so that the previous 
convergent Pn-i/q-n-i was "hit" at the end of a segment), Im means that it combines with the 
previous partial quotient in a block R^'^'^JL (so that pn-i/q-n-i was "missed"), and Ic means 
that it combines with the previous partial quotient in a block with a C-symbol (a "corner"). 
However it does not specify which sequence of symbols actually occur as legal expansions. We 
study this next. 

5. Additive Continued Fractions to Vertical Cutting Sequences 

In this section, we show how to construct the Minkowski geodesic continued fraction expan- 
sion of 9 given the additive continued fraction expansion of 9. Let 9 have the ordinary continued 
fraction expansion 

9 = [oo; 01,02, 03, • • •], 

in which oq = or — 1. In view of the results of section 4, it suffices to determine for each 
On = 1 whether or not it has label l/j, 1^ or Ig in the partition of the MGCF expansion given 



in Theorem 4.1. 



Theorem 5.1. Given 9 = [oq; oi, 02, 03, • • •] with — ^ < ^ < ^, suppose o„+i = 1. Set an = 
[0,o„,o„„i, . . . ,ai] = ^tlzl^ and [in = [on+i, an+2, • • so that an G [0,1] and (3n S [1,2]. In 



In-l 

terms of the linear fractional transformation 



N(z) := 



1 2 

2 1 



(z) = , (5.1) 

^ ^ 2z + l ^ ^ 



the Minkowski geodesic continued fraction expansion of 9 has 



if (3n > N(q„) , 

an+i = { Ic if Pn = N(an) , (5.2) 
if (3n < N(an) . 
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Remark. The matrix N acting as a linear fractional transformation maps [0,1] to [1,2] 

-1 2 



is not integral. The theorem 



while reversing orientation. Its inverse N ^ '^2 1 

could also be formulated in terms of the linear fractional transformation — N which maps [0, 1] 
to [—2,-1] and is an involution, i.e. — N(— N(z)) = z. In particular, since N sends M"^ into 
itself, 

/3 < N(a) ^ a > -N(-/?) . (5.3) 

The symbol Ic can occur only when Q is rational, because N(a„) is always rational, while /3n is 
rational if and only if the expansion of Q terminates. 

Proof of Theorem |5.1| . Note that if ai = 1, it always becomes 1^, and we have a„ = 
and thus N(q„) = 2 > /3„ as required. The following discussion assumes that On+i = Ij with 
n > 1, so that we are not in the segment Bq. 

By the argument in the proof of Theorem 4J, whether a„+i becomes 1/^, Ic, or 1^ depends 
on which of ( ^I^ ) and ( [4.6D holds with equality for a larger value of t\ we have \m if it is (^^), 
and Ic if both hold with equality for the same t. Since a^+i = 1, we have 



Pn+l = Pn+Pn-1- 



The condition for 1^ is thus 



{Pn - qnOf - {Vn+l " Q'n+l^ ^ {Pn-1 - Qn- lOf " {Pn " Q'n^')^ 



9n+l In 



or, equivalently. 



yn^*)^ - {Pn+1 - Qn+lOf 



> 



In-l 



In Qn+1 



(Pn^l - Qn-lOy - {Pn " qnOY In-l " In 

and the condition for Ic is equality. Putting everything in terms of Pn-i, Pn, Qn-i, and Qn gives 



(5.4) 



(5.5) 



{Pn-l - qn-lO) + 2 (p„_i - qn-l9) {pn - QnO) Qn-l + "^Qn-lQn 



{Pn-l - Qn-lO) - {Pn 



Now let 



> 



9n-l In 



(5.6) 



a. 



Pn 



In-l 



(5.7) 
(5.8) 

These quantities have ordinary continued-fraction expansions a„ = [0, a„, a„_i, . . . ,ai], and 
Pn = [an+i>«n+2> ' ' ']; See Venkov section 2.4]. Clearly < < 1 and /5„ > 1, and, since 
a„+i = 1, Pn < 2. Now (|5.6|) becomes 



Qn 

Pn-l - 



/3n - 2/3n al + 2a, 
/32-1 ^ al-1 



(5.9) 
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or, by subtracting 1 from each side, 



2/3„ + l 2a„ + l 

-^3^ > ^J^- (5-10) 



02-1 a2 - 1 

Let C = (2a„ + 1)/(q^ — 1). Then ( ^.lOf ) holds if and only if — is between the roots of 

Cx'^ -2X-1-C = {). (5.11 



The sum of the roots of ( |5.11| ) is 2/C, and the larger root is a^, so the other root is 

2a^ - 2 _ -a„ - 2 ^ -N(aO- (5.12) 

This gives our desired condition; we have a„+i = 1^ if /3„ < N(a n). Since — N is an involution 
and is increasing on [—2, —1] and on [0, 1], we can also write this condition as an > — N(— /3n). 



Theorem |5.l| suffices to classify all 9 containing the symbol Ic, i.e. all vertical geodesies that 



hit a corner of a translate of a fundamental domain. 

Corollary 5.1. Let — ^ < 6 < ^ and suppose that the MGCF expansion of contains a symbol 
Ic- Then 6 is rational and has ordinary continued fraction expansion of the form 

[0,ai,a2,,...,an,lc,6i,---,6m]> «/0<6'<i. ,^ 

[-1, 1, ai - l,a2, ,an, Ich, ■■■ ,bni], if-^<9<0. 

in which the additive continued fraction expansion of a*^ := [l,bi, . . . ,bm] is computable from 
the additive continued fraction expansion of On = [0, a„, . . . ,ai] by a finite automaton, and 
vice-versa. Furthermore there is an absolute constant cq such that 

ACF-length «) < cq (ACF-length (an)). (5.14) 



Proof. Raney |43] proves that given any fixed linear fractional transformation M(z) := 

with integer coefficients a, b, c, d and ad — bc^ 0, the ACF expansion of M(^) can be computed 
from the ACF expansion of 9 using a finite automaton. Furthermore the conversion process 
inflates the ACF-length by at most a multiplicative constant (depending on M). Apply this to 
Pn = N(q„). The specific automaton for N(z) = is given in Raney pp. 274-275], and 
the constant in this case is 3. I 

These results give a bound on the computational complexity of computing the MGCF 
expansion from the additive continued fraction expansion. 

Theorem 5.2. The Minkowski geodesic continued fraction expansion of 9 can be computed 
from the additive continued fraction expansion in quadratic time using linear space. That is, 
there are absolute constants ci, C2, C3 such that for any 9 with —^<d<^, the first £ symbols of 
the MGCF expansion of 9 can be computed using the first cii symbols of the additive continued 
fraction expansion of 9 using at most c^^ time steps and c^i space locations. 



27 



Proof. Theorem 4.1 shows that the main problem is to resolve whether a given symbol 
a„+i = 1 which appears as the ith symbol in the ordinary continued fraction expansion of 9 
is to be Im, l/i! or Ic in the MGCF expansion. For this, we use Theorem At worst, we 
must look all the way back to the beginning of the MGCF expansion. Corollary shows that 
the ACF expansion for N(a„) is of length at most cqI, and thus comparing it to (3n requires 
looking at no more than coi + 1 symbols, and this uses 0{l) time and 0{l) space. Since there 
may be 0{i) different I's to be resolved, the total time required is 0{i'^). Note that we test 
the inequalities ( p.2\ ) by comparing symbol sequences for N(a„) with the initial part of that 
for f3n- This algorithm can easily be implemented on a random access machine with the given 
time and space bounds. It can also be implemented on a one-tape Turing machine with the 
same space bound and a time bound polynomial in i. We omit details. For these two standard 
computational machine models, see Aho, Hopcroft and Ullman [^. I 

There are examples which do require time steps; for example, 

e = (V3-l)/2 = [0,2^1], 

has this property. For this 6, it is necessary to backtrack all the way to the first symbol to 
determine that each 1 is Ih, because the sequence [0, (2, l/i)-', 2, Ic, (2, 1^)-', 4] is a Ic-sequence 
for each j > 1- 

In the next section, we will show that the MGCF expansion cannot be computed from 
the additive continued fraction expansion using a finite automaton, or even using a pushdown 



automaton with one stack, as defined in Hopcroft and Ullman |21|. 



6. Vertical Geodesies: Forbidden Blocks 

In this section, we characterize the allowable cutting sequences of 11^. 

Definition 6.1. A finite word W in the symbol set {L, R, J} is a forbidden block of 11^ if it 
occurs in no cutting sequence of 11^; otherwise, it is an admissible block. It is an excluded initial 
block of if does not occur as an initial segment of any cutting sequence of 11^; otherwise, it 
is an included initial block. 



It is easy to see that excluded initial blocks alone can be used to characterize any subset 
of the one-sided shift on {L,R,J}. We show that 11^ has the stronger property that it is 
determined by its set of forbidden blocks, as follows. 

Theorem 6.1. A finite word W in the symbols L, R, J is an excluded initial block in if and 
only if at least one of the blocks LW and RW is a forbidden block ofH^. 

Proof. We prove the contrapositive. First, we show that for every included initial block 
W, both R"^W and L"W are admissible blocks of 11^ for all n > 1. Geometrically, this encodes 
the fact that the geodesic {oo,6) is a limit of geodesies {9',0) as \6'\ — > oo, where 9' — oo is 
associated to R^W and 9 +oo is associated to L"W. The block W corresponds to a finite 
initial segment of the cutting sequence 9 + it for some irrational 9, say for to ^ ^ ^ oo and it 



has first symbol J. By Lemma 2.2, this geodesic cannot hit a corner of any translate of hence 



there is some positive e such that it is at distance at least e from any corner for to < t < 9. Pick 
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and thus produces the sequence W. Now take a matrix M 



G SL{2,Z), so that 
'he PSL{2, Z)-action 



p/q a large negative rational number, and observe that the geodesic {0',0) contains the word 
R"W in its one-sided infinite cutting sequence, in which it enters the domain ^ just after the R". 
Here we require that the radius r = 9 — p/q is at least n + 2 so that it produces the sequence R , 
and satisfies r — \/r'^ — 1 < e so that the geodesic is within e of (oo, 6*) over the range to < t < 1 

q' p' 
q -p 

M(7) = (oo, 9') with 9' = M(0), and note that 9' is necessarily irrational, 
doesn't affect cutting sequences, so the cutting sequence of M(7) still contains the word R W. 
There is a unique choice of q' and p' such that — ^ < ^' < ^. This vertical geodesic has the 
same cutting sequence as 7. Thus R W and hkewise L W are admissible blocks of 11^ for all 
n > 1. 

Second, we prove that if W = Si • • • S„ is a word in L, R, J for which RW and LW are both 
admissible blocks of 11^, there is a vertical geodesic whose cutting sequence has W as an initial 
segment. To show this, note that the first symbol in W is necessarily J, since LR and RL are 
forbidden blocks. Let 71 and 72 be irrational vertical geodesies whose cutting sequences contain 
the words RW and LW, respectively. Translate each of them by the appropriate elements of 
PSL{2,Z) to geodesies 7J = {pi/qi,9[) and 72 = (^2/921^2) such a way that the translated 
geodesies enter the fundamental domain at the symbol immediately preceding W in their 
cutting sequences. Since the cutting sequence of 7^ has the letter R as it enters it is oriented 
in the direction of increasing real part and has pi/qi < —1/2, while 72 is oriented in the 
direction of decreasing real part and has P2/92 > 1/2. The first symbol in W is J, hence each 
geodesic exits through its bottom edge at a point with real part between —1/2 and 1/2. If 
7^ and 72 do not intersect, that is if 9[ < 92, then pick any irrational 9 between 9[ and 92, say 
9 = 9[, and let 7 = (00, 9). If j[ and 72 do intersect, let 9 be the real part of their intersection 
and let 7 = {00, 9). In either case, 7 passes through J" since < 9[ < 9 < 9'2 < ^. Thus 
the intersection must have real part 9 with —^<9<^. In either case, we show that 7 has 
initial word W in its cutting sequence. After leaving J^, both j[ and 72 pass through the same 
sequence of translated fundamental domains J^i, ■ ■ ■ ,Tn, entering along the same edge of each, 
corresponding to the symbols in W. The vertical geodesic 7 must then enter Tj on the same 
edge at a point Zj between the points zij and Z2^j where 7^ and 72 hit it, because the domains 
Tj are hyperbolically convex; see Figure 6.1. Also, it cannot hit any other fundamental domain 
between and Zj, because every point in that interval on 7 is on a geodesic between points 
of 7^ and 72 which are in T . Thus the initial cutting sequence of 7 is W. Finally, if the finite 
endpoint ^ of 7 is rational, then since 7 hits no corners up to hitting the edge W„, it is a 
Euclidean distance at least e away from every corner on the W-edges. For every irrational 9' 
with \9' — 9\ < e, the irrational vertical geodesic 7' = {9' + it : t > 0} has the same initial 
segment W. I 

There are two main types of forbidden blocks (including initial blocks). Some blocks are 
edge- forbidden; they cannot occur in a cutting sequence because they correspond to a geodesic 
hitting two edges which are part of the same hyperbolic line. Others are whole-forbidden because 
a geodesic can hit any pair of edges in the sequence, but no geodesic can hit the whole sequence. 

Series |4^, Theorem 3.1] shows that any minimal edge- forbidden block for any geodesic flow 
must correspond to geodesic hitting a sequence of domains which are all adjacent to the same 
boundary line, crossing that line as the first and last step. This allows us to find the edge- 
forbidden blocks from Figure 2.3; without loss of generality, we can let the boundary line be the 
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ii 7 72 




P2 g Pi 

92 qi 

Figure 6.1: Intersecting geodesies: 7 is trapped between 7^ and 72 within J^j. 
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unit circle, with the forbidden block starting inside the unit circle, going outside, and coming 
back inside. Since the unit circle has three segments, there are nine such blocks: JJ, LR, RL, 
□□, RJRJ, JUL, JRJR, LJLLJL, RJRRJR. These also lead to excluded initial blocks. Since 
the edges R and L meet at infinity, any forbidden block beginning with either one corresponds to 



an excluded initial block. These blocks are exactly those which are forbidden by Theorem 4.2 



We next obtain from Theorem 5.1 a characterization of the the whole-forbidden blocks and 
whole-excluded initial blocks in 11^. These conditions involve a critical symbol Ih or 1^ which 
is mislabeled in a cutting sequence expansion. 

Definition 6.2. A finite sequence [di, . . . ,dn, 1, 6i, . . . ,bm] of positive integers is ambiguous if 
it does not suffice to determine the MGCF-label on the critical symbol 1. More precisely, if 

6o := [0,d„,...,(ii], and 5i := [0, d„, . . . 1], 
(3o ■■= [l,6i, . . . ,6m], and /5i := [1, 6i, . . . ,6^ + 1] , 

then [di, . . . ,dn, l,bi, . . . is ambiguous if and only if 

int[/3o,/3i] nint[N(5o),N(5i)] / 0, (6.1) 

Here int[5o)'^i] denotes the closed interval determined by 6o and 6i with any ordering of the 
endpoints, i.e. 6i < 6o may occur. 

Note that the interval int[N((5o)) N((5i)] contains all real numbers whose continued-fraction 
expansions begin with the same partial quotients as Sq. 

Definition 6.3. A central sequence [di, . . . , dn, l,bi, . . . ,bn] is an ambiguous sequence such 
that p/q = [0 or — l,di, . . . ,dn, l,bi, . . . ,bm] is a Ic sequence. We say that p/q is the rational 
number associated to this central sequence. 

We use the following definition to designate initial words of using the initial symbol 
di = oo to indicate that there is no preceding term in the cutting sequence but that the 
geodesic comes from oo. 

Definition 6.4. An initial sequence [oo, d2, ■ ■ ■ , dn, 1, 6i, . . . , bm] denotes an initial word, which 
if — ^ < ^ < has d2 = 1, ds = ai — 1, and di = ai-2 for i > 4, while ii < d < ^, di = Oj_i 
for ah i > 2. 



We define ambiguous initial sequences as in Definition |6.2| . If we set 

5o = 5i = [0,dn, ... ,d2] = [0,dn, . . . ,d2,oo], (6.2) 
then the initial sequence is ambiguous if 

int[/3o, n int[N(5o), N((^i)] / 0, (6.3) 
where int[N((5o), N(5i)] is the single point N(5o)- 
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Theorem 6.2 (Characterization of 11^) The set 11^ of cutting sequences of irrational 9 in 
— ^ < ^ < ^ uses the alphabet {R, L, J}. It consists of all sequences that factorize in segments 
B0B1B2... as in Theorem \4-^ with the additional property that no sequence of consecutive 
segments corresponds to a Minkowski geodesic continued fraction expansion 

[di, . . . , d„, 61, . . . , bm] 

with di = 00 allowed, with associated 5o, 61, (3o, Pi, such that: 

(i) Both of [^2, . . . , dn, 1, 61, ... , bm\ and [tii, ^2, . . . , dn, 1, 61, ... , b„i-i\ are ambiguous, 
(a) Either 

U = Ih and int[/3o,/3i] < int[N(,5o), N(5i)], (6.4) 

or 

U = Im and int[/3o,/?i] > int[N(,5o), N(5i)] . (6.5) 



Proof. Immediate from Theorem using Theorems and I 

Theorem 6.2 characterizes a large set of forbidden blocks which are given by conditions (i) 
and (ii) when di is finite. It also gives an additional set of excluded initial blocks when di = 00. 
These are sufficient to determine 11^ because they identify each 1 as Ih or 1^ as appropriate. 
However, Theorem |6.2| does not give a complete set of minimal forbidden blocks, because there 
are extra forbidden blocks in which several symbols di = 1 and bi = 1 are replaced by Ih and 
Im, which if left as indeterminate I's would not be forbidden. 

The complete set of minimal forbidden blocks of 11^ seems harder to characterize. Consider 
a given block of symbols L, R, J; this can be parsed by Theorem into a block of complete 
segments [ai, . . . , am] in which each symbol is either 1^, Im, or an integer at least 2, together 
with some conditions on the adjacent incomplete segments; for example, an incomplete segment 
L can be part of a complete segment encoding either Om+i > 4, or Om+i = 3, am+2 = Im- Then 
[oi , . . . , Om] is a forbidden block if and only if a certain finite set of linear fractional conditions 
on two real numbers a, (3 of the form 

a < — (6.6) 

cf3 + d ^ ^ 

with ad — be = ±1 are inconsistent. There is one such inequality for each symbol 1^ or 1^ 
in the block, which encodes the condition that [a,ai, . . . ,am, (3] produces the correct symbol 
Ih or Im- There may also be one or two inequalities on a alone (or on /3 alone) if there is an 
incomplete segment at that end, and one more condition (|6.6| ) if the possible encoding of that 
segment includes another symbol Ih or Im- The total number of inequalities is linear in m. 

Tables 6.1, 6.2, and 6.3 below list some central sequences, ambiguous sequences and for- 
bidden blocks, respectively; these were obtained by applying the transformation N to simple 
6i- Table 6.3 illustrates the computation of forbidden sequences from a central sequence; any 
change to the terms in a central sequence which increases either a„ (and thus decreases N(a„)) 
or f3n forces the 1 to be Ih, and conversely for Im- 



Theorem 6^ implies that 11^ is complicated in the sense that its set of minimal forbidden 



blocks is very large. 
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[2,1c, 4] 
[1,1,1c, 4] 


[2,1c, 3,1] 
[1,1,1c, 3,1] 


[3, lc,2,2] 
[1,2,1c, 2, 2] 


[3,1c, 2, 1,1] 
[1,2, lc,2,l,l] 


[4, lc,2] 
[1,3,1c, 2] 


[4, lc,l,l] 
[1,3,1c, 1,1] 


[3i + l,lc,l,i] 

[i,3j, ic,i,i] 


[3i + l,lc,l,i-l,l] 
[l,3j, lc,l,i-l,l] 


[3j + 2,lc,l,j,3] 
[l,3j + l,lc,l,j,3] 


[3j + 2,lc,l,j,2,l] 
[l,3j + l,lc,l,j,2,l] 


[2, 2,1c, 3] 
[1,1, 2,1c, 3] 


[2,2, lc,2,l] 
[1,1, 2,1c, 2,1] 


[2, 3,1c, 2, 5] 
[1,1, 3,1c, 2, 5] 


[2,3, lc,2,4, 1] 
[1,1, 3,1c, 2, 4,1] 


[3, 2,1c, 3, 4] 
[1,2, 2,1c, 3, 4] 


[3,2, Ic, 3, 3,1] 
[1,2, 2,1c, 3, 3, 1] 


[4, 3,1c, 2, 3] 

[1,3,3, lc,2,3] 


[4,3, lc,2,2,l] 

[1,3, 3,1c, 2, 2,1] 


[5, 2,1c, 3, 2] 
[1,4, 2,1c, 3, 2] 


[5, 2,1c, 3, 1,1] 
[1,4, 2,1c, 3, 1,1] 


[j,l,lc,3j + l] 

[l,j-l,Llc,3i + l] 


[j,l,lc,3j, 1] 
[l,j-l,l,lc,3j, 1] 


[3,j,l,lc,3i + 2] 
[l,2,j,l,lc,3j + 2] 


[3,i,l,lc,3i + l,l] 
[l,2,j,l,lc,3i + l,l] 



Table 6.1: Some central sequences, including all with at most five terms. 



[•••,J,1,1, 3i + l,any,...] 


[...,any,3j + l,l, l,i, ...] 


[•••,J,l,l,3i + 2,any,...] 


[...,any,3j + 2,1, l,j, ...] 


[•••,i,l,l,3j + 3,any,...] 


[...,any,3j + 3,l, l,j, ...] 


[...,1,2,1,2,1,...] 




[...,2,2,1,3,>4,...] 


[...,>4,3,1,2,2,...] 


[...,3,2,1,3,2,...] 


[...,2,3,1,2,3,...] 


[ . . . , 3, 2, 3, 3, . . . ] 


[...,3,3,1,2,3,...] 


[...,4,2,1,3,2,...] 


[...,2,3,1,2,4,...] 


[...,>5,2, 1,3,1,...] 


[...,1,3,1,2,>5,...] 



Table 6.2: Non-central ambiguous sequences which go two terms forward and two terms back 
from the underlined 1. 
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[> 3, > 


3] 




[> 3, Ih, 2, 


1] 




[1,1m, 1] 






[> 1,2,1„,1] 






[1, Im, 2, 


> 


2] 


[> 1,2,1,„,2, 


> 


2] 



Table 6.3: Forbidden sequences obtained from the central sequence [1, 2, Ic, 2, 2]; the reverses 
of these sequences are also forbidden. 

Theorem 6.3. The number n{k) of minimal forbidden blocks o/H^ of length at most k grows 
exponentially in k; in fact 

liminfn(A:)^/'' > 2^/^2_ (g 7) 

fc— »oo 

Proof. We will show that each central sequence associated to a rational p/q ^ 1/2 yields 
two minimal forbidden blocks. The forbidden blocks are produced by adding one symbol to 
each end of the central sequence, and by replacing the central Ic with a three-symbol word. 
All these forbidden blocks are distinct. Assuming these facts are proved, consider the central 
sequences 

[di,d2,..., dn, Ic, bi,..., km]. (6.8) 

in which each dj = 1 or 2, and in which [bi, . . . ,bm] is determined from [di, . . . ,dn]- The 
number of symbols in the cutting sequence encoding of [di, . . . ,dn] is at most 3n, and by 
Corollary the number of symbols in the cutting sequence encoding of . . . , bm] is at most 
9n. In obtaining forbidden blocks, the Ic term is encoded by three symbols, and one symbol is 
added at each end, hence all resulting forbidden blocks contain at most 12n + 5 symbols. We 
conclude that there are at least 2""^^ minimal forbidden blocks of length at most 12n + 5, which 
proves ( |6.7D . 

We now construct the minimal forbidden blocks. We use the result of Appendix A, which 
shows that any vertical geodesic for 9 = p/q with — | < ^ < ^ hits at most one corner of a 
fundamental domain. This fact implies that each of the central sequences ( |6.8D contains only 
one symbol Ic, and thus each p/q comes from at most one central sequence. Associated to the 
central sequence is a word W1SW2 in which S = Ci or C2 is a corner symbol. For the corner 
symbol, there are two choices of three symbols in {R, L, J}, which replace the Ic by 1/j or 1^; 
for Ci, the choices are JRJ or LJL. We consider the eight words obtained from W1SW2 by 
replacing S by either choice of a three-symbol block, and by adding a prefix symbol and a suffix 
symbol, each of which may be either L or R. The continued fraction for p/q can be recovered 
from any one of these eight blocks, and since it has only a single Ic, the critical Ic in the central 
sequence is also uniquely determined. 

We claim that six of these eight blocks are admissible blocks for 11^ and the other two are 
forbidden blocks. The vertical geodesic 7 = {00, p/q) corresponding to 

p 

- = [0 or - l,di, ...,dm, Ich, • • • ,M 
Q 

can be approximated by geodesies in the symbol topology in six different ways. Two of these 
consist of approximating geodesies which do not cross 7 at all but approach it from the left 
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and right, respectively. The other four consist of approximating geodesies which cross 7, either 
crossing above or below the corner that 7 hits at the Ic, and initially approaching 7 either from 
the left or from the right. For sufficiently good approximations, these produce the admissible 
blocks. (Here we again use the fact that the geodesic 7 hits exactly one corner, which implies 
that all geodesies sufficiently close to 7 have the same sequences Wi and W2.) The other two 



blocks are forbidden by condition (ii) of Theorem 6.2. They encode cutting sequences for a 
geodesic that would have to approach 7 from the left, pass to the right of the corner, and then 
return to the left of 7, or vice versa; such a geodesic would have to cross 7 twice, which is 
impossible. 

Let C(j and C;, denote the two possible three-symbol encodings of S. Of the four blocks 

RWiCaWzR , RWiCaWsL , LWiCaWaR and LWiCaWaL , 

exactly three are admissible and one is forbidden. Every sub-block of the forbidden block 
appears in one of the three admissible blocks, hence the forbidden block is minimal. The 
same argument applies to the other four blocks containing C;,, and produces a second minimal 
forbidden block. I 

Theorem 6.4. There does not exist a finite automaton which, when given the additive contin- 
ued fraction expansion of an irrational number 9 with Q < < \ as its input sequence, computes 
the Minkowski geodesic continued fraction expansion of 6, or, equivalently, the cutting sequence 
expansion of the geodesic {00, 9). 

Proof. Any finite automaton that would compute the Minkowski geodesic expansion of 9 
must output the n-th term of this expansion after seeing at most a bounded number of symbols 
following the n-th symbol of the ACF expansion of 9. However the forbidden block criteria of 
Theorem |6.2| show that it is sometimes necessary to see an arbitrarily large string of symbols 
after the n-th symbol, to decide if Ih or 1^ should be used. For example, for each j > 1, the 
sequence 

[24^-+2,le,3,(8,4y] (6.9) 

is a central sequence. Adding some later terms decreases /34j+2j and thus changes the Ic to 
Im, while decreasing the final 4 in (8,4)-' to a 3 and then adding some further terms increases 
/?4j+2 and thus changes the Ic to l/^. A finite automaton cannot look ahead through the 14j -|-6 
terms which are necessary, since j can be any integer. In fact, even a pushdown automaton 
(with one stack) cannot correctly compute all such cutting sequences, because it must look 
14j -|- 6 steps ahead, but also 12j -|- 6 steps back to distinguish [oq; 2'^-', 3, (8, 4)-^ , 8, 3, . . .] 
from [ao; 24^+2,^,3,(8, 4)^8, 3,...]. ■ 



7. Two-sided Cutting Sequences: Structure of Sjr 

We now use the information on vertical cutting sequences H^ to characterize the two-sided 
cutting sequences T,jr. 

Theorem 7.1. The cutting sequence shift T,jr for the fundamental domain T of PSL{2,'Ij) is 
the closed subshift whose forbidden blocks coincide with the set of forbidden blocks o/H^. 



35 



Proof. This result follows from the fact that the set of images of the set of irrational vertical 
geodesies under PSL{2,Z) is dense in the space of all geodesies of S^/PSL{2,Z). 

Let and 5£ and 5s denote the complete sets of forbidden blocks of 11^, S^r, and Sjr, 
respectively. Since Sjf is the closure of S^, we have 5^ = S^. 

We first show that 

5° c 5°. (7.1) 

For this it suffices to show that every admissible word W in 11^ is an admissible word in S^. 
Suppose that the word W appears in the cutting sequence of the irrational vertical geodesic 
(oo, 9) = {9 + it : t > 0}. This geodesic is a limit of geodesies {(j)i,9) where ^ oo through a 
sequence of values such that {(pi, 9) hits no corner of an image of J-. The word W = Si . . . 
corresponds to a specific set of edges of translated fundamental domains {gjJ^ : < j < r} 
with Qj € 5L(2,Z). For all sufficiently large the geodesic {(pi, 9) passes through the same 
sequence of fundamental domains {gjJ' : < j < r}, hitting the same sequence of edges in the 
same order. Thus W occurs in the two-sided cutting sequence of {(pi, 9), so it is an admissible 
word of S^. 

The reverse inclusion 

50 C Si (7.2) 

is proved similarly. Let W be an admissible word in some general position geodesic 7 = (</>, 9) . 
There is a sequence of translated fundamental domains {J-'j : < j < r} which (p passes 
through that corresponds to W. Since 7 hits no corners, we can choose a rational number 
p/q sufficiently close to (p such that the geodesic 7' = {p/q,9) passes through the same set of 
fundamental domains {J-'j '■ < j < r}, hitting the same sequence of edges in the same order, 
hence its cutting sequence (which is only one-sided infinite) contains the word W. Now apply 
to 7' the transformation M in PSL{2,'L) which takes p/q to cxd, and 9 to some 9' in {—\, \)- 
For this choice, the cutting sequence of M(7') is in H^, and ( [7.2| ) follows. I 

Theorem 7.2. Every element in Sjf is a cutting sequence for a unique oriented geodesic on 
9)/ PSL{2,'L) which hits T , except j or the two sequences R°° or L°°. Every oriented geodesic 7 
on 9)/ PSL{2,'L) has at least one and at most finitely many shift- equivalence classes of cutting 
sequences in Sjf. 7/7 is not periodic then it has at most eight shift- equivalence classes of cutting 
sequences in Sjf. 

To establish this result, we first prove a preliminary lemma. 

Lemma 7.1. Let jj = {9'p 9j) for j = 1,2, . . . be a sequence of general position geodesies that 
intersect T which have cutting sequences {S^^ : i G Z} such that the symbol Sq^ corresponds to 

the geodesic jj entering the fundamental domain T . If the cutting sequences {S^^} converge in 
the symbol topology as j ^ co to a limit sequence {Sj} then the endpoints 9'- and 9j converge 
to unequal limiting values 

9' = lim 9'. and 9 = lim 9j. (7.3) 

in MU {—00, cxd}. The geodesies 7j converge to a limiting geodesic 7 = {9', 9) if at least one of 
9,9' is finite. The exceptional cases (—00, cxd) and (00,— cxd) correspond to the limiting symbol 
sequences and L°° , respectively. 
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Proof. Let S+ = {Si : i > 0} and S_ = {Si : i < 0} denote the positive and negative part 
of the hmit sequence, respectively. If S+ is R°°, then since any geodesic which enters J- after 
So and has S^'''' = R for 1 < i < n must have 6j > n + 1/2, we conclude that linij^oo 9j = co 
in this case. Similarly, if S+ is L°°, then Ivcuj^ao (^j = — oo; if S_ is R°°, then limj^oo ~ 
if S_ is L°°, then limj^oo 9j = —oo. In particular, this shows that the two-sided limit sequence 
R°° corresponds to (—00,00), and L°° corresponds to (00,— 00). 

We next observe that since RL and LR are forbidden blocks, if S+ is not R or L it must 
contain a symbol J. Let the initial segment of S+ up to the first such symbol be R"J (resp. 
L^J). Any geodesic 7^ which matches these symbols necessarily has 9j < n + 1 since it crosses 
the semicircle with endpoints n — 1 and n + 1 (resp. 9j > —n — 1) and is positively oriented 
(resp. negatively oriented), hence 6j > —1/2 (resp. 9j < 1/2). In either case, if a limit ( ^ ) 
exists it must be finite. Similar reasoning applies to Sj to show that limj^oo d'j is infinite if and 
only if S_ does not contain a symbol J. 

We now suppose that S+ contains a symbol J, and claim that limj^oo Gj exists and is finite. 
The argument above shows that {6j} is bounded, so to prove this claim it suffices to show that 
if {9j} is not a Cauchy sequence then the one-sided sequences {Sp'' : i > 0} cannot converge in 
the symbol topology. If it is not a Cauchy sequence, there is some e > such that for any A^, we 
have j,k > N with \9j — 9k\ > e. The general-position geodesies jj and 7^ enter and thus have 
radius at least y/3/2. Consequently, any point z = x + iy on the geodesic 7j with \x — 9j\ < 1/4 
and < y < e' for e' < e^/36 will actually have \x — 6j\ < e/3, and thus \x — 9k\ > 2e/3; 
similarly, points on 7^ with < y < e' and \x — 9^1 < 1/4 must have \x — 9k\ > 2e/3. This 
implies an upper bound on the number of fundamental domains which 7^ and 7^ both hit 
somewhere inside the box \Re{z) — {9j + 9ii.)/2\ < 2e, Im{z) < e. We may assume e < 1/2, hence 
any such fundamental domain has a cusp at a finite rational point p/q. The rosette R{p/q) 
of all translated fundamental domains that touch p/q has Euclidean diameter less than 4/g^. 
Thus any translated fundamental domain J^' which is hit by both jj and jk inside the box 
must have two points whose real parts differ by at least e/3, hence q < yl2/e. In addition, p/q 
must lie within Euclidean distance 4/g-^ of both 9j and 9k, hence there are only finitely many 
such rational p/q. Finally, there is a finite bound depending only on q and e on the number 
of fundamental domains in each rosette R{p/q) which have Euclidean diameter exceeding e/3. 
Together this yields a finite upper bound on the number of common fundamental domains 
J^' between 7^ and 7^ in the box. Outside the box, after exiting the geodesies 7j and 7^ 
traverse the strip e' < Im{z) < 1, and thus can hit at most N'^ translated fundamental domains 
J^' in that strip, because there are only finitely many geodesies with Euclidean radius at least 
e' which intersect the region — ^ < Re{z) < i. We conclude that the positive cutting sequences 
7j and 7^ cannot agree on all of their first A^ -|- A^ -|- n -|- 1 symbols. Thus the sequences 
S^^ = {S,^^^ : j > 0} cannot converge in the symbol topology, which proves the claim. 

A similar argument applies to the limit sequence S_ if it contains the symbol J, to prove 
that lim,-^oo = d' exists and is finite. 

The geodesic 7 = (0', 9) is encoded as a limit of the sequence of geodesies 7^. If all of the ji 
for i > N intersect a particular domain J^' , then since the endpoints O'^ and 9i converge to 9' and 
9, and is closed, 7 also intersects J^, at least on the boundary. If no 7^ for i > N intersect a 
domain J^' , then 7 does not intersect the interior of JF'. Thus, if we consider 7 to hit a domain 
if either it passes through the interior, or it intersects the boundary and all 7^ for sufficiently 
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large i intersect the interior, we can define a cutting sequence as the set of edges which separate 
these domains, and this sequence is the limit in the symbol topology of the cutting sequences 
of 7i. ■ 



Proof of Theorem 7.2, Since Sjp is a compact set, every oriented geodesic 7 = {0' ,9) 
which hits the interior of J- has at least one cutting sequence in C(7) by taking the limit point 
of cutting sequences from a family 7^ = {OpOj) with 9j — > 9' and 0j — > as j — > 00. 

It remains to show that each oriented geodesic in f) which hits J- has only a finite number of 
shift-equivalence classes of cutting sequences in [C(7)]. The simplest cases are general position 
geodesies, for which C{'~f) is a single cutting sequence, and [C(7)] consists of a single shift- 
equivalence class. 

We determine how many shift-equivalence classes of cutting sequences are possible for a 
limiting geodesic. First, consider the limiting geodesies with exactly one rational endpoint 
(considering 00 as rational). These correspond to vertical geodesies under the PSL{2, Z) action, 
so by Lemma such geodesies cannot hit any internal corners, and thus have two possible 
encodings; either can be obtained by using approximating geodesies which approach the rational 
endpoint from either side. Thus [C(7)] contains two shift-equivalence classes. 

Next, consider those geodesies which have two rational endpoints. Without loss of gener- 
ality, we can move one endpoint to 00 and so have a vertical geodesic 7 = {00, p/q). Such a 
geodesic hits only finitely many corners of fundamental domains. If it hits n corners at finite 
values of t, then C{^) contains exactly 2n + A cutting sequences. Two of them come from 
geodesies approximating 7 from either side without crossing it, and the other 2n + 2 result from 
approximating geodesies which cross 7 between the kth. and {k + l)st corners with < < n, 
either from left to right or from right to left. (Here the 0th corner is 00, and the {n + l)st 
corner is p/q.) In Appendix A, we prove that n < 1 for all rational p/q, except those p/q = \ 
(mod 1), which have n = 2. Thus [C(7)] contains at most eight shift-equivalence classes. 

It remains to bound the number of shift-equivalence classes of cutting sequences for geodesies 
which have two irrational endpoints. If such a geodesic hits no corners, then it has only one 
cutting sequence in C(7). If it hits exactly one corner, then it has exactly two cutting sequences, 
which are obtained using geodesies which approach it while staying on opposite sides of the 
corner. The difficult ease occurs with geodesies that hit at least two corners. We show that all 
such geodesies hit infinitely many corners and are periodic. In this case, C{'y) will be an infinite 
set. To show periodicity, we observe that a corner in the upper half-plane is the intersection of 
two circles with centers at rational points on the x-axis and rational radii, so its x-coordinate 
is rational and its y-coordinate is the square root of a rational number. Thus, if the circle with 
radius r and center (xo,0) passes through two such points (xi,yi) and (x2,y2), then it satisfies 
the equations 

{xi-xof + yl = r^, 
{x2-xof + y2 = r^. 

Since yf and 2/2 rational, so is r^; also, equating the left sides gives a linear equation 
for xq with rational coefficients. Thus the circle intersects the x-axis in two algebraically 
conjugate real quadratic surds. Pell's equation allows us to find M G 5'L(2,Z) which preserves 
one endpoint (and thus its conjugate). Applying this transformation to our geodesic re-scales 
it while preserving the orientation; hence the geodesic necessarily is a periodic geodesic on 
Sj/PSL{2,'Z). The set C{'y) is then infinite because we can choose a set of approximating 
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geodesies which cross between any pair of consecutive corners of 7, and the resulting hmit 
cutting sequences are all distinct. 

We complete this case by showing that if a periodic geodesic 7 on Sj/PSL{2, Z) hits exactly 
n corners of translates of fundamental domains in its period, then there are exactly 2n + 2 
shift-equivalence classes in [(^(7)]. Label the corners which are hit in one period ci,...,Cn- 
The approximating geodesies 7, can only be close in the symbol topology if they all cross 7 
between the same pair of corners, or if none cross 7 at all. If they cross, then we can use the 
periodicity of 7's cutting sequence to shift the crossing point between ci and c^+i. With n 
possible crossing regions, and crossing possible either from inside to outside or vice versa, there 
are 2n encodings; we get two more from approximating geodesies which are completely inside 
or outside 7, for a total of exactly 2n + 2 shift-equivalence classes. I 

Theorem gives no uniform upper bound on the number of shift-equivalence classes of 
cutting sequences that correspond to periodic geodesies. We formulate the question whether a 
uniform upper bound exists as open problem (3) in the concluding section. There are periodic 
geodesies which correspond to 10 cutting sequence shift-equivalence classes because they hit four 
corners in one period, such as (— \/l3, V^) and (— -v/133, Vl33). Could this be the maximal 
number that occurs? 



Theorem 7.1 implies that the shift Sjp is a relatively complicated set. Indeed Theorem |6.3| 
showed that the set of minimal forbidden blocks is very large. We now show that Sjp is not 
a sofic shift. A sofic shift is any shift that is a factor of a shift of finite type (see Marcus and 
Lind IQ].) Alternatively, it is the set of possible bi- infinite walks on an edge-labeled finite 
graph. 

Theorem 7.3. The cutting sequence shift Ejf for the fundamental domain T of PSL{2, Z) is 
not a sofic shift. 

Proof. If S is a shift, the follower set i<s(W) of a one-sided word W = (. . . ,S„2)S_i) is 
the set of all one-sided words = (So, Si, . . .) such that WW"*" is an element of S. Sofic shifts 
are characterized by the property that the totality of different follower sets {i*i;(W) : all W} is 
finite; see Marcus and Lind p4| . Theorem 3.2.10]. 

We show that Sjf has infinitely many follower sets. The sequences [3, 2^-'"'^^, Ic, 3, (8, 4)-' , 10] 
and [4,2''J+2,lc,(8,4y,13] are central. Thus, in [. . . , 3, 24^+2, U, 3, (8, 4)^ 13, . . .], the U is Ih 
if j > k but Im if j < k, regardless of the symbols on either side. Thus the follower sets of 
one-sided words which end [3, 2^^'^'^] are different for all j. I 

8. The Shift Ejr Determines JF up to Isometry 

Our object is to prove the following result. 

Theorem 8.1. Let T he a finitely generated discrete subgroup of PSL{2,M) that acts properly 
dis continuously on Sj, and which has a polygonal fundamental domain V which is hyperbol- 
ically convex. Suppose that the cutting sequence shift S-p is isomorphic to '^j^^psL(2,z) ^ 
permutation of symbol alphabets. Then there is an element g G PSL{2,M) such that T = gV 
and 

PSL{2,R) = gVg-' . 
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Proof. Since the cutting sequence contains only three symbols, the polygon must be a 
triangle. 

Since a geodesic can hit the L edge an unlimited number of consecutive times, the L and 
L edges of the triangle must intersect at an angle of zero. The L edge cannot be the 
same as the L edge because a geodesic can hit it twice consecutively, nor can it be the J edge 
because no geodesic can hit that edge twice consecutively; thus it must be the R edge. Since a 
geodesic cannot hit the J edge twice consecutively, the J edge must equal the J edge. Thus 
the generator at the J edge is an involution, and to have determinant 1, it must be an inversion. 

Let the angle between J and L be a, and the angle between J and R be /3. If 2a + /? < tt, 
then a geodesic approaching the JR-corner could hit edges J, L, J, L in sequence; however, this 
cutting sequence is not possible for our fundamental domain. Likewise, 2/3 + a > vr. Also, for 
any hyperbolic triangle, the sum of the angles is less than vr, so a + /? < vr. 

The conditions for a given polygon to be a fundamental domain are given in Maskit p5| , 
section 2]. In our case, it is necessary that either n{a + P) = 2tt for some n, or na = n and 
mP = TT for some m and n. The only cases consistent with the conditions on a and /? are 
a = P = it/3, the desired fundamental domain; and a = it/3, /3 = tt/2, which is the half of our 
fundamental domain with x > 0. But that domain is only a fundamental domain for a group 
including the reflection in the line x = 0, which has determinant —1. (That domain also has 
different symbolic dynamics; a geodesic can hit L and R consecutively.) I 



9. Open Problems 

(1) . We have shown in one special case that the polygon V can be recovered from the data 
S-p up to isometry. The example comes from a Riemann surface of genus 0. Does the result 
persist for higher-genus Riemann surfaces? Can any S-p be explicitly determined for a Riemann 
surface of genus at least one? 

(2) . Since Sjf determines up to isometry, in principle it determines vol{S^/PSL{2,Z)). 
Can vol{S^/PSL{2,Z)) = ^ be easily computed directly from Sjf? 

(3) . Is there a universal upper bound on the number of shift equivalence classes of cutting 
sequences corresponding to any periodic geodesic on S)/ PSL{2,'L)1 Equivalently, is there a 
universal upper bound on the number of times that a periodic geodesic can hit a corner of the 
fundamental domain during a single period? 

(4) . The zeta function C,t.{z) of a shift S is defined by 

OO h. 



Ce(^) =exp(J^iVfc — 



A,•=l 



in which counts the number of periodic words in S of period k. Is there a simple formula 
for the (dynamical) zeta function of Sjf? What is the topological entropy of Ejf? 

(5). Cutting sequence shifts S-p can be constructed in higher-dimensional cases along the 



lines considered in |gO[, if one restricts to a suitable subclass of geodesies, called "flat" geodesies 
in [^]. Can any such Sp be determined explicitly? 

Acknowledgment. We are indebted to L. Flatto and M. Sheingorn for helpful comments 
and references. 
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A. A Bound on the Number of Corners on a Vertical Geodesic 



Lemma A.l. For a rational 0, the vertical geodesic ^ = {6 + it : t > {f\ has at most one value 
oft such that 9 + it is a corner of a PSL{2,'Z) -translate of T , unless Q = \ (mod 1), in which 
case it has exactly two such values, which are t = V3/2 and Vs/Q. 

Proof. Since the corner — 1/2 + \/—3/2 of J- is obtained from the corner 1/2 + \/—3/2 by 
the transformation z — > z — 1, every corner is a PS'L(2, Z)-translate of 1/2 + -v/~3/2. Let the 

a b 



element of SL{2, 12) be 



d 



with ad — be = 1 . We then have 



a b 
c d 



1 V- 

2^~ 



a|i + V 



+ b 



d 



(4ac + 2ad + 26c + ibd) + {2ad - 26c) 
4(c2 + cd + d2) 

(2ac + ad + 6c + 2bd) + 

2(c2 +cd + d2) • 



(A.l) 



Let iV = 2ac + ad + 6c + 26d and D = c^ + cd + d"^ , so that the real part of ( |A.1D is N/2D. 
We win show that N/ 2D is either in lowest terms or can be reduced to lowest terms by dividing 
both and D by 3. (The factor 3 can occur; for example, take a = 2, 6 = 1, c = 1, d = 1.) 
Since ad — bc = 1, ad + 6c is odd and thus N is not divisible by 2. Now substitute a = (l + 6c)/d 
in ( [A.iP ; this gives 

N _ 2^c + 2bd+^d + bc 
2D ~ 2D 

_ 2c + 2bc^ + 2bd'^ + d + 2bcd 
2d(c2 + cd + d2) 

6(c2 + cd + d2) + 2c + d 



2d(c2 + cd + d2) 

Since (c, d) = 1, we have c^ + cd + d'^ = 1 (mod 2), and also (d, c2 + cd + d2) = (d, c2) 
follows that 



(A.2) 



1. It 



{N,D) 



{dN, D) = {bD + 2c + d,D) = (2c + D,D) = (2c + d, c2 + cd + d2) 
= (2c + d,c^ + cd + d^ - d(2c + d)) = (2c + d, c(c - d)). 



Now, (c, 2c + d) = (c, d) = 1, and (2c + d, c — d) = (3c, c — d) < (3, c — d)(c, c — d) = 1 or 3. 
Therefore (A^, Z?) = 1 or 3. 

Suppose the geodesic is 6 + it with 6 = v/w. There are only two possible values of ( [A.l| ) 
on this geodesic; N/2D can only equal v/w in lowest terms if w = 2D ox w = 2D/2). The 
imaginary part of ( |A.1| ) is either \/—3/uj or \/—3/3w. However, these two points are separated 
by a hyperbolic distance of In 3, the length of the finite side of the domain J^; therefore, they 
can both be PSL{2, Z)-translates of corners only if they are connected by a PSL{2, Z)-translate 
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of an edge. This edge is the geodesic 9 + it itself, and the only vertical edges of translates gj^ 
occur when 9 = (n + 1/2) for n G Z. 

For 9 = n+1/2, there are two SL{2, Z)-translates of corners. If we take a = 1, b = n, c = 0, 
d = 1, we get the corner (n + 1/2) + -y—3/2; if we take a = n + 1, 6 = n, c = 1, c/ = 1, we get 
the corner (n + 1/2) + ■ 
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